A  SEARCH  FOR  COSMIC  AXIONS 


By 

CHRISTIAN  A.  HAGMANN 


A  DISSERTATION  PRESENTED  TO  THE  GRADUATE  SCHOOL 

OF  THE  UNIVERSITY  OF  FLORIDA  IN  PARTIAL  FULFILLMENT 

OF  THE  REQUIREMENTS  FOR  THE  DEGREE  OF 

DOCTOR  OF  PHILOSOPHY 


UNIVERSITY  OF  FLORIDA 
1990 


UNIVERSITY  OF  FLORIDA  LIBRARIES 


ACKNOWLEDGEMENTS 

It  is  my  great  pleasure  to  thank  my  advisors  Profs.  N.  Sullivan  and  D. 
Tanner  for  their  guidance,  valuable  advice  and  general  support  during  the  time 
of  my  research.  I  am  grateful  to  Prof.  P.  Sikivie  for  his  continued  interest  in 
my  research  project  and  many  stimulating  and  helpful  comments.  It  was  a  great 
honor  for  me  to  participate  in  the  Florida  Axion  Project  and  become  a  member 
of  a  new  and  exciting  branch  of  physics. 

My  thanks  go  also  to  Profs.  G.  Bosman,  H.  Campins,  J.  Ipser,  and  M.  Meisel 
for  their  interest  and  patience  while  serving  on  my  committee.  The  experiment 
would  not  have  been  possible  without  the  support  by  the  cryogenic  group  under  A. 
Hingerty  and  G.  Labbe,  the  machine  shop  under  B.  Fowler  and  the  electronics  shop 
under  L.  Phelps.  I  am  grateful  to  C.  Porter  for  providing  the  3-D  plotting  program. 
I  would  also  like  to  mention  my  fellow  graduate  students  and  the  staff  members 
of  the  Physics  Department  for  creating  a  pleasurable  working  environment. 


TABLE  OF  CONTENTS 

Page 

ACKNOWLEDGEMENTS ii 

ABSTRACT    v 

CHAPTERS 

1  INTRODUCTION 1 

2  ORIGIN  AND  PROPERTIES  OF  THE  AXION      5 

The  Axion  and  the  Strong  CP  Problem 5 

Axions  in  Astrophysics  and  Cosmology 14 

Axions  as  a  Dark  Matter  Candidate      25 

Detection  of  Galactic  Axions 29 

3  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 37 

Magnet  and  Cryogenic  System 37 

The  Cavities 40 

The  Microwave  Amplifiers 57 

Room  Temperature  Data  Acquisition  Hardware 66 

4  DATA  ACQUISITIION  AND  ANALYSIS     73 

The  Search  Process 73 

The  Data 75 

Limit  on  the  Electromagnetic  Coupling  of  the  Axion 79 

5  SUMMARY  AND  CONCLUSION 83 

APPENDICES 86 

A     CAVITY  DEVELOPMENT 86 

iii 


Page 

B     SOFTWARE  LISTINGS 93 

REFERENCES 120 

BIOGRAPHICAL  SKETCH 125 


Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the     , 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 


A  SEARCH  FOR  COSMIC  AXIONS 
By 

Christian  A.  Hagmann 

August  1990 

Chairman:  Neil  S.  Sullivan 
Cochairman:  David  B.  Tanner 
Major  Department:  Physics 

A  search  experiment  has  been  carried  out  to  look  for  relic  cosmic  axions 
trapped  in  the  halo  of  our  Galaxy  and  limits  on  their  abundance  are  presented 
for  the  axion  mass  range  of  5.48  x  10~^  eV  <  rua  <  6.77  x  10~^eV  and  7.47  x 
10-^  eV  <ma<  7.60  x  10"^  eV. 

The  detector  consisted  of  a  high-Q  microwave  cavity  immersed  in  liquid 
helium  and  permeated  by  a  strong  magnetic  field.  A  low  noise  microwave  receiver 
with  a  cryogenic  HEMT  amplifier  as  its  first  stage  was  coupled  to  the  cavity.  The 
noise  power  within  a  cavity  bandwidth  was  spectrum  analyzed  and  searched  for 
narrow  peaks  resulting  from  axion  to  photon  conversions  via  the  Primakoff  effect. 
The  resonance  frequency  was  changed  by  displacing  a  dielectric  rod  inside  the 
cavity. 


CHAPTER  1 
INTRODUCTION 


The  axion  is  a  hypothetical  particle  which  was  proposed  to  explain  why  the 
strong  interactions  conserve  the  symmetries  of  P  and  CP.  The  Lagrangian  of 
QCD  contains  a  symmetry  violating  term 

^^  =  ^Tr(G^.Gn  (1-1) 

Id  tt 

where  the  G^i/'s  are  the  gluon  fields,  g  is  the  QCD  coupling  constant  and  6  e 
(0,27r),  is  a  free  parameter  of  the  theory.  To  be  consistent  with  experiments,  6 
has  to  be  very  small  (9  <^  10"^) ,  which  seems  unnatural. 

Peccei  and  Quinn  [1]  presented  an  elegant  solution  to  this  problem  by  in- 
troducing a  special  spontaneously  broken  global  [/(I)  symmetry  into  the  QCD 
Lagrangian.  This  broken  symmetry  has  the  effect  of  making  6  a  dynamical  vari- 
able, which  adjusts  itself  to  the  CP  conserving  value  ^  =  0.  As  pointed  out  by 
Weinberg  [2]  and  WiJczek  [3],  there  mtist  exist  a  Goldstone  boson  associated  with 
the  spontaneously  broken  U{1)pq  symmetry;  it  was  named  the  axion. 

Massless  at  the  classical  level,  the  a:xion  acquires  a  small  mass  due  to  quan- 
tum effects.    The  mass  was  originally  thought  to  lie  in  the  keV  range,  making 

it  accessible  to  accelerator  searches  and  beam  dump  experiments.  However,  this 
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particular  search  was  unsuccessful,  and  a  new  kind  of  'invisible'  a^ion  was  pro- 
posed [4,  5].  This  new  axion  has  a  much  smaller  mass  and  is  very  weakly  coupled 
to  ordinary  matter.  Mass  and  couplings  of  the  axion  are  suppressed  by  the  PQ 
symmetry  breaking  scale  vp(^,  with  the  mass  given  by 

„„  ,  lO-^V  (i5^!°^)  (1.2) 

V         VpQ         J 

and  with  all  couplings  proportional  to  rria- 

Invisible  in  most  terrestrial  experiments,  axions  can  have  important  impli- 
cations for  astrophysics  and  cosmology.  Thermally  produced  axions  can  freely 
stream  out  from  the  interior  of  stars  and  contribute  to  their  radiative  energy  loss. 
In  order  not  to  cool  stars  excessively,  the  axion  couplings  need  to  be  sufficiently 
small,  and  this  constraint  places  an  upper  limit  on  the  axion  mass.  An  even  better 
limit  comes  from  the  recent  observation  of  supernova  SN1987a.  Axion  emission 
can  dominate  nevitrino  emission  for  a  small  axion  mass.  Yet  the  observed  neutrino 
pulse  can  account  for  all  of  the  energy  released  during  the  collapse.  This  can  be 
used  to  place  an  upper  boimd  on  the  axion  mass  of  ma  ^  10~    eV. 

Cosmological  considerations  provide  a  lower  limit  on  ma-  When  the  temper- 
ature in  the  early  universe  falls  below  the  PQ  scale,  U  (l)pQ  gets  broken  sponta- 
neously and  topological  defects  in  form  of  cosmic  strings  will  appear.  This  string 
network  will  oscillate,  causing  the  emission  of  axions,  which  contribute  to  the 


present  energy  density  of  the  universe  and  woiold  close  it  for  rria  ^  10~^*  eV. 
If  an  inflationary  period  occurs  after  the  PQ  phase  transition,  the  strings  will 
become  diluted  and  no  axions  are  produced  from  strings.  However  in  that  case, 
axions  are  produced  due  to  the  initial  misahgnment  of  the  vacuum  angle  6. 

For  TUa  ~  10~^  eV,  cosmic  axions  from  string  radiation  or  initial  vacuum 
misalignment  constitute  a  significant  fraction  of  the  energy  density  and  would  be 
all  or  part  of  the  dark  matter  which  is  known  to  exist  in  our  universe.  In  particular 
our  own  galaxy  has  a  dark  halo  of  unknown  composition  which  could  conceivably 
be  made  out  of  axions. 

The  experiment  described  herein  has  made  an  attempt  to  detect  these  galac- 
tic axions  for  an  axion  mass  near  10~^  eV.  As  first  proposed  by  Sikivie  [6],  the 
detection  scheme  was  based  upon  the  electromagnetic  coupling  of  the  axion  field 
a,  given  by 

-^077  =  ^077  O  E  •  B  ,  (1-3) 

with  ga-ff  being  proportional  to  iria-  In  a  static  magnetic  background  field  Bq,  this 
coupling  stimulates  the  conversion  of  axions  into  photons.  Bq  provides  a  virtual 
photon  and  a  real  photon  carries  off  all  the  axion  energy.  Since  galactic  axions 
are  nonrelativistic,  the  photon  energy,  uj  =  ma(l  +  O(10~^)),  will  have  only  a 
small  spread.  The  frequency  of  these  photons  lies  in  the  microwave  region,  and  the 


decay  rate  of  axions  is  resonantly  enhanced  in  a  high-Q  microwave  cavity  tuned 
to  the  axion  energy. 

The  cavity  is  kept  at  hquid  heUum  temperature  to  reduce  its  thermal  emis- 
sion and  is  coupled  to  a  very  sensitive  microwave  receiver.  The  signature  of  an 
axion  conversion  signal  would  be  an  increase  of  the  power  inside  of  the  cavity 
above  the  noise  floor.  The  cavity  has  to  be  tunable  in  order  to  search  a  wide  range 
of  possible  axion  frequencies.  Tuning  is  accomplished  by  moving  dielectric  and 
(or)  metal  rods  inside  of  the  cavity,  yielding  tuning  ranges  of  about  10-20  %.  To 
date,  two  cavities  have  been  constructed,  which  enabled  us  to  cover  the  frequency 

range 

1.32GHz  </<  1.63GHz     and   1.80GHz  </<  1.83 GHz  (1.4) 

corresponding  to  axion  mass  ranges  of 

5.48  X  10-^eV  <  m^  <  6.77   x  IQ-^eV     and 

(1.5) 

7.47  X  IQ-^eV  <ma<  7.60  x  10" W  . 
No  positive  signal  coming  from  axion  conversion  has  been  detected  so  far, 
and  this  result  enables  us  to  place  an  upper  limit  on  the  axion  coupling  for  the 
above  mass  range,  assuming  that  the  galactic  halo  is  axionic. 


CHAPTER  2 
ORIGIN  AND  PROPERTIES  OF  THE  AXION 

The  Axion  and  the  Strong  CP  Problem 

Quantum  chix)modynamics  (QCD)  ,  the  gauge  theory  of  strong  interactions, 
sucessfuUy  describes  many  aspects  of  quark  physics  such  as  scahng  and  asymptotic 
freedom.  The  Lagrangian  is  usually  written  as 

k 

/here 


1  ^ 

L  QCD  =  "o  Tr  G^uG^'''  +  J2^k  (n'^D/z  -  r>n)^k  (2-1) 


G^iy  =  d^A„  -  ^l,A^^  -  ig  [A^j,,Ai,] 

Dpq^  =  {d^  -  ig  A^)  q;t  ^2.2) 

8 
A,  =  Y.AI\^I2 

a=l 
and  where  the  G^v's  are  the  gluon  fields,  the  Aa's  are  the  Gell-Mann  matrices 
and  k  runs  over  all  quark  flavors  :  u,d,s,c,t,b  .  The  classical  QCD  Lagrangian 
possrases  a  high  degree  of  symmetry  a,nd  in  particular  conserves  P,  CP  and  T. 

Due  to  the  non-Abehan  nature  of  the  gluon  fields,  there  exist  vacuum  states 
of  topologically  distinct  finite  energy  field  configurations  [7]  labeled  by  a  winding 
number  n.  Upon  quantization,  tunneling  solutions  (instantons)  appear,  and  the 
true  vacuum  state  becomes  a  linear  combination  of  the    \n )  states.   The  ground 


state  has  the  form  of  a  Bloch  wave  and  is  given  by 

l^)  =  E^"''"'l")-  (2.3) 

n 

The  '6-%'acvitim'  modifies  the  Lagrangian  that  enters  the  path  integral  formulation 
of  QCD  by  adding  a  new  term  ; 

Leff  =  L  QCD  +  ^Tr  (G^.G'-)  (2.4) 

where  G^'^  =  T^e^^^^P  Gj^p  and  is  the  dual  of  G^*^.  The  symmetry  properties  of 
GG  are  the  same  as  of  FF  in  electrodynamics  and  are  odd  under  P,  CP  and  T. 
The  coefficient  6  is  an  arbitrary  constant  with  modulus  27r,  and  the  effective 
Lagrangian  violates  P,CP  and  T,  provided  ^  ^  0,7r.  The  best  known  physical 
manifestation  of  the  ^-vacuum  would  be  a  nonvanisliing  electric  dipole  moment  of 
the  neutron,  dn-  The  current  upper  bound  of  dn  is  [8] 

\dn\  ^  10"^^  e  cm  '         (2.5) 

which  constrains  6  to  be  very  small,  namely 

1^1  ^  10-^  .  (2.6) 

The  strong  CP  problem  is  a  question  of  how  to  make  6  naturally  small. 


The  situation  becomes  even  more  mysterious  if  the  electroweak  interactions 
are  included.  The  spontaneous  symmetry  breaking  of  the  electroweak  SU{2)  ® 
U{1)  gauge  group  is  responsible  for  giving  mass  to  the  quarks.  The  mass  matrix 
in  the  Lagraagian  is  in  general  off-diagonal  and  non-hermitian,  i.e. 

L  mass  =  -qm  Mij  qij  -  qii  i^^)ij  ^Rj  ■  (2-7) 

One  can  diagonalize  the  mass  matrix  Mij  with  two  special  unitary  transformations 
acting  on  the  left  and  right  quark  fields  respectively, 

qu  -^  Uiij  qij 

(2.8) 

qm  -*  Ufiij  qnj , 

leaving  the  mass  matrix  in  the  form  M,j  =  |mjj|  <5jj  e'".    The  common  phase 
factor  a  can  be  removed  by  a  chiral   17(1)^  transformation  with 


qii  -  ^-^'"Z'  qii 


qm  -^  e'"^^  ^Ri 


(2.9) 


bringing  the  mass  term  into  its  final  form 


L  mass  =  -  X]  "^Jt  qifc  qik  •  (2-10) 


The  latter  transformation  affects  the  QCD  angle  6  through  the  well  known 
Adler-Bell-Jackiw  (ABJ)  anomaly  [9]  of  the    U{1)a  axial  current: 

^5  =  ^^kl^  75qifc 

k  (2.11) 

2 
dJ^  =  2nt  -^Tr  {GnvG^^)  +  mass  term  contributions    . 

The  effect  of  performing  a  chiral  rotation   a/2  is  to  change   6  according  to 

^  -^  ^  +  «_;•  a  =  6"  +  arg  det  M  =  'S  .  (2.12) 

'B  is  the  effective   CP  violating  parameter  in  the  QCD  Lagrangian,  and  it  replaces 

6  in  Eqs.  (2.4)  and  (2.6). 

The  most  convincing  solution  of  the  strong  CP  problem  was  put  forward 
by  Peccei  and  Quinn  [1].  They  postulated  the  existance  of  a  new  spontaneously 
broken  global  U{\)pQ  symmetry,  thereby  making  'O  a  dynamical  variable  which 
adjusts  itself  naturally  to  zero.  The  U{\)pq  chiral  symmetry  is  spontaneously 
broken  at  a  scale  vpn,  and  there  exists  a  Goldstone  boson  called  the  axion  [2]. 
The  axion  field  transforms  as  the  phase  of  a  complex  scalar  field: 

a  —*  a  -\-  vpqSa  (2-13) 

where  6a  is  the  infinitesimal  parameter  of  the  U{\)pq  transformation.    Exact 
Goldstone  bosons  have  only  derivative  coupUngs,  but  here  it  is  postulated  that 


the  U{\)pQ  symmetry  suffers  from  the  same  ABJ  anomaly  as  the  ordinary  QCD 
U{1)a  symmetry.  This  anomaly  results  in  a  potential  for  the  axion  and  gives  it  a 
mass,  making  it  a  pseudo-Nambu-Goldstone  boson.  The  modified  Lagrangian  is 


^     ^\    3 


Leff=L  QCD  +  (<?  + OrKTr  (G^^Cn  +  -d^a  d^'a  +  L  (5^a,  i>)  (2.14) 

'  ^  ^  VpQ  IblT^  ^ 

where  ^  is  a  model  dependent  constant  and  L(d^a,\f)  specifies  the  couplings  to 
fermions.  The  axion  potential  is  minimized  for  [1] 

{a)  =  -~QvpQJi  (2.15) 

which  is  the  vacuum  expectation  value  of  a.  The  physical  axion  field  corresponds 
to  excitations  about  the  minimum: 

(^physical  =  a-{a)  .  (2.16) 

The  first  axion  model  embedded  the  U{l)pQ  quasi-symmetry  in  the  standard 
electroYi?eak  theory.  The  standard  Higgs  doubkt  #  was  replaced  by  two  separate 
doublets  [10] 

y/2     '  \0j  '         %/2 


1  m/l\  ^  1  m/0\  ,„,^, 

^,=         .,e^[\         $2  =  ^7=^2-     J  (2.17) 


where 

x='^         v  =  yjvj  +  v^  =  {V2Gp)-2  -  250  GeV  =  vpg  (2.18) 
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and  where  Gp  is  the  Fermi  constant.  The  vacuum  expectation  values  v\  and  U2 
break  down  the  symmetry  SU{2)yY  ®  UO-)y  ®  UO-)pQ  '^'^'^^  ^0-)em  giving  rise  to 
four  Goldstone  bosons.  Three  of  them  are  absorbed  through  their  couphngs  to  the 
gauge  fields  generating  mass  for  the  vector  bosons  W''^  ^W~ ,  Z  .  The  remaining 
Goldstone  boson  is  the  axion. 

Current  algebra  techniques  are  used  to  calculate  the  mass  and  the  lifetime 
of  the  axion.  The  result  is  [2,  10,  11] 

ma  =  rriT: TT  [^  +  -) 


vpQ    2  X    1  +  z  ^2.19) 

~  12n/(x  +  -)keV 
■'  X 

where  ttit^  and  f-^  are  the  pion  mass  and  decay  constant  and  z  =  mxi/m^  ~  0.56. 
For  rria  <  2me  the  dominant  decay  channel  is  a  -^  27  with  the  lifetime 

Ta:^  { sec.  2.20 

\     TUa     / 

This  kind  of  axion  has  been  seardied  for  unsucessfully  in  accelerator  experi- 
ments. One  important  par2Lmeter  is  the  branching  ratio  for  kaon  decay  into  axions 
which  is  calculated  to  be  [12] 

BR(A'+ ^7r+  +  a)>  1.3x  10"^  (2.21) 


11 


while  the  experimental  upper  hmit  is  3.8  x  10~°  [13].  Another  negative  result 
comes  from  heavy  quarkonium  decays  where  the  gq-^  a  +  7  process  has  not  been 
observed. 

Not  long  after  these  negative  results,  new  axion  models  were  proposed  mak- 
ing the  axion  very  weakly  interacting  and  impossible  to  see  in  accelerator  exper- 
iments. Two  models  of  the  'invisible'  axion  are  usually  considered  by  particle 
theorists.  The  KSVZ  or  hadronic  axion  [4]  leaves  the  standard  electroweak  model 
intact  and  introduces  a  heavy  exotic  quaxk  Q  and  a  complex  scalar  field  a,  both 
of -whidi  carry  Peccei-Quinn  charges  and  are  singlets  under  SU{2)\y  <S)  U(1)y  ■ 
Under  the  PQ  symmetry,  they  change  according  to 

(2.22) 

,     -2iQ„ 

(T  — >  e         a 

while  leptoas  and  ordinary  quarks  are  singlets.  The  Yukawa  coupling  between  Q 
and  a  is 

Ly  =^-fQL<^QR-  f*QR^*  Ql  (2.23) 

and  the  modified  Higgs  potential  has  the  form 

V($,  a)  =  -^l  $t$  _  ^la*a  +  A$($t$)2  +  \^(a*af  +  A^.^^^^  ^V       (2.24) 

where  $  is  the  Higgs  doublet  from  the  standard  model.  Ly  and  V'($,<t)  are  left 
invariant  under  U{l)pQ. 
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In  addition  a  discrete  symmetry 

R:        Ql^-Ql,  Qr-^  +Qr  ,  cr-.-a  (2.25) 

is  required  to  avoid  a  bare  mass  term  mQQ  .  The  a  field  develops  an  arbitrary 
large  vacuum  expectation  value  (a)  =  vpQ  above  the  electroweak  scale,  and  the 
phase  of  a  becomes  the  axion  field  a.  The  axion  couples  to  ordinary  quarks 
and  photons  indirectly  through  Q  and  gluon  loops.  There  is  also  a  higher  order 
coupling  to  electrons  via  radiative  corrections.  The  axion  mass  is  given  by 


m 


KSVZ_^_,    _/^^«.-,n-6.vfl^!!^£X)  (2.26) 

V       VPO       / 


~6x  10~W( 


1  +  ^  VpQ  V        VpQ 

while  the  photon  coupling  is 

'  ga-r^^aFt"'Ft'''^gajjaE-B,  (2.27) 


with 


_       a       /E      2  A  +  z\ 
~  2-Kvpo  ViV  ~  3  1  +  2/' 


KSVZ 


— 0-GeV-'(^-il±f)(,^) 


(2.28) 


where  a  ~  ^  is  the  fine  structure  constant  and  E/N  is  the  ratio  of  electromag- 
netic and  color  anomaly  of  U(1)pq. 
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The  DFSZ  axion  [5]  model  has  two  Higgs  doublets  as  in  the  original  model. 
In  addition  there  exists  a  scalar  field  a  analogous  to  the  one  in  the  KSVZ  model 
with  VEV  vpQ  >  V.  Both  o  and  ordinary  fermions  carry  PQ  charges.  The  a 
field  couples  indirectly  to  quarks  and  leptons  via  a  quartic  interaction  between  a 
and  $i,$2?  with  the  latter  fields  having  Yukawa  interactions  with  fermions.  The 
relevant  interaction  between  the  Higgs  fields  is  given  by 

l/($i,  #2,  o)  =  \  ($7*2<^^)  +  h.c.  (2.29) 

which  fixes  the  PQ  transformation  properties  of  a  relative  to  $i,$2  ^^^d  the 
fermions.  In  this  scenario,  the  mass  and  couplings  are  similar  to  the  KSVZ  axion 

with 

DFSZ  \A    _       /:r  /„  on^ 

•'   1  +  2  VpQ 

and 

nFSZ_„      Q     /8      2A  +  z\ 

^«^T'  ^vpq\Z      Z\^z)  ^2  31) 

=  1.38xlO--GeV-(^^). 
The  PQ  symmetry  breaking  scale  vpQ  for  both  KSVZ  and  DFSZ  axions  is  con- 
strained from  above  only  by  the  GUT  scale    mpianck  —  ^0^^  GeV,  corresponding 
to  a  mass  ma  —  10~      eV. 
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Axions  in  Astrophysics  and  Cosmology 

Astrophysical  considerations  provide  powerful  constraints  on  invisible  axion 
parameters.  Axion  emission  contributes  to  the  energy  loss  of  stars  and  affects 
their  evolution.  The  energy  transport  out  of  stars  in  the  form  of  axions  is  maxi- 
mized for  axion  couplings  yielding  a  free  mean  path  of  roughly  one  stellar  radius. 
Increasing  the  coupling  strength  leads  to  a.xion  trapping  inside  an  axion  sphere 
whereas  decreasing  the  coupling  allows  free  streaming  but  reduces  the  number  of 
axions  produced.  In  order  not  to  exceedingly  shorten  observed  stellar  lifetimes, 
cooling  through  axions  must  be  a  sufficiently  weak  process  and  thus  yields  limits 
on  axion  couplings. 

In  main  sequence  stars  and  red  giants  the  dominant  processes  for  axion 
coohng  are  Compton-like  scattering  [  e~  +  7  ->  e~  +  a  ]  and  bremsstrahlung  [  e~  + 
Z  (e~)  — »  e~+  Z  (e~)  +a  ],  whereas  for  KSVZ  axions  with  no  tree  level  couphng  to 
leptons,  the  most  important  loss  mechanism  is  the  Primakoff  effect  [f+Z  (e~)  — > 
a  +  Z{e~)  ].  The  associated  Feynman  diagrams  are  shown  in  Fig.  2.1 .  AppHed 
to  red  giants,  these  processes  lead  to  an  upper  bound  on  the  axion  mass  in  the 
KSVZ  model  of  [14] 


15 


/ 


/ 

/a 


/ 


J- 


b) 


^ 


/ 


N 


/ 


^ 


N 


/ 
/ 


7T 


7 


N 

C) 


N 


Z.e 


Fig.  2.1:  Feynman  diagrams  of  axion  production  mechanisms  in  stars  with  a) 
PrimakofF  effect,  b)  Compton-hke  scattering,  c)  nucleon-nucleon  bremsstrahlung 
and  d)  electron-nucleon  and  electron-electron  bremsstrahlung. 
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In  the  DFSZ  model,  axion  cooling  will  prevent  red  giants  from  entering  the  helium 
burning  phase  unless  [15] 

^DFSZ  <;^Q-2gy  (2.33) 

For  an  axion  mass  of  m^  >  O  (10— 300keV) ,  the  previous  arguments  are  no  longer 
valid  because  thermal  axion  production  is  suppressed  by  the  Boltzmann  factor 
e~T^.  However,  this  mass  region  is  already  excluded  by  laboratory  experiments 
[13]. 

The  most  recent  astrophysical  limit  comes  from  supernova  SN  1987a.  Dur- 
ing core  collapse,  the  enormous  gravitational  binding  energy  of  ~  10^^  erg-  was 
released  and  radiated  away.  According  to  the  standard  picture,  this  energy  was 
carried  away  in  the  form  of  a  neutrino  pulse.  Due  to  the  high  density  and  temper- 
ature of  the  core,  neutrinos  are  trapped  inside  a  neutrino  sphere  and  are  only  ra- 
diated from  the  surface.  Axions  with  sufficiently  small  mass  and  couplings  stream 
freely  and  couJd  greatly  accelerate  the  cooling  process  if  produced  abundantly.  The 
dominant  production  process  for  axions  in  supernovae  is  nucleon  bremsstrahlung 
[]V  +  iV  — >  iV-l-iV  +  a]  (see  Fig.  2.1)  with  a  cross  section  proportional  to  rua- 
For  increasing  axion  mass,  trapping  occurs  and  the  cooling  efficiency  decreases. 
Therefore  one  ends  up  with  an  axion  mass  range  for  which  axion  emission  domi- 
nates or  is  comparable  with  neutrino  emission.  The  neutrino  pulse  from  SN  1987A 
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[16]  was  consistent  with  a  purely  neutrino-cooled  core,  which  yields  the  excluded 

mass  region  [17] 

10-3eV<  ma  ^  2 eV,  (2.34) 

which  is  valid  for  the  KSVZ  and  DFSZ  model.  The  combination  of  astrophysical 
bounds  allows  a  mass  range  for  DFSZ  axions  of 

m?FSZ<io-3eV  (2.35) 

and  two  windows  for  the  KSVZ  axion  of 

m^SVZ  ^  10-3  eV     ^^^     2eV  <  m^^Z  ^  g^y  .  (2.36) 

Another  important  bound  on  the  axion  mass  comes  from  cosmological  considera- 
tions. In  the  early  universe,  axions  may  have  been  in  thermal  equihbrium  with  the 
primordial  plasma.  A  rough  criterion  for  a  particle  to  be  in  equiUbrium  is  T  >  /T 
where  F  is  the  interaction  rate  with  other  particle  species  and  H  is  the  universal 
expansion  rate  or  Hubble  constant.  Because  the  interaction  rates  between  axions 
and  other  particles  scale  hke  the  axion  mass,  only  axions  heavier  than  ~  10~  eV 
will  be  in  thermal  equilibrium  at  some  time.  These  thermally  produced  axions  will 
eventually  fall  out  of  equilibrium  as  the  the  universe  expands  and  their  number  per 
comoving  volume  will  remain  constant  except  for  the  decay  channel  [a -^  2f].  The 
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contribution  that  these  axions  add  to  the  present  energy  density  of  the  universe 
is  given  by  [18] 

^t/ierma/^^^  10"^  (^a/eV,)  (2.37) 

where  fi  =  pipe  and  100 /ikmsec"^  Mpc~^  is  the  present  value  of  the  Hubble 
constant.  From  Eq.  (2.37)  it  is  clear  that  thermal  axions  could  make  up  a  signifi- 
cant fraction  of  the  present  energy  density  only  for  mass  ranges  which  are  already 
excluded  by  stellar  evolution. 

A  possibly  much  more  important  contribution  [19,  20]  to  the  energy  density 
comes  from  the  misalignment  of  the  initial  B  value  with  the  CP  conserving  value 
of  ^  =  0.  When  the  temperature  of  the  early  universe  falls  below  the  PQ  breaking 
scale  vpQ,  the  Higgs  field  a  takes  on  the  form 


a  =  vpn  e 


MhpQ  ■  (2.38) 


where  the  phase  assumes  a  random  value  in  the  range  0  <  a  <  27r  vpQ  and  the 
axion  field  will  in  general  be  a  function  of  position  a  =  a  (x) .  However  if  a  period 
of  inflation  occurs  during  or  after  the  PQ  phase  transition,  a  small  patch  of  space 
will  grow  exponentially  during  that  epoch  and  render  a  constant  in  the  observable 
part  of  universe.  As  the  temperature  of  the  universe  reaches  T  ~  1  GeV,  instanton 
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effects  become  important  and  the  axion  acquires  a  small  mass  given  by  [20] 

ma(T)~  0.1  rr7a(T  =  0)(^^)^■^     T  »  Aqcd  (2-39) 

where  Aqcd  —  200  MeV  and  is  the  QCD  scale,  below  which  quarks  are  confined 
in  hadrons.  Below  Aqcd  the  axion  mass  will  assume  its  value  given  by  Eqs.  (2.26) 
and  (2.30).  A  finite  axion  mass  provides  a  potential  for  a  of  the  form  (see  Fig.  2.2) 

Via)  =  ^ma{Tfa'^  +  0{a'^)  (2.40) 

which  drives  the  axion  field  to  the  CP  conserving  value  a  —  0.  Actually,  in  some 
axion  models,  V{a)  can  have  N  degenerate  minima  (^  =  A''  in  Eq.  (2.14))  which 
can  be  modeled  as 

na)  =  ">I(T)(^)'(i-co.(^))  ,2.«) 

and  can  lead  to  a  domain  wall  problem  [21].  In  the  DFSZ  model,  N  =  6. 

When  the  axion  mass  ma{T)  becomes  larger  than  the  Hubble  factor  H{T), 
the  axion  field  begins  to  oscillate  at  a  frequency  determined  by  ma{T).  The 
Lagrangian  governing  the  axion  field  dynamics  is 

L=^{d^af-V{a),  (2.42) 

yielding  the  equation  of  motion 

^  (a)  +  3  Hit)  ^  (a)  +  mlit)  {a)  =  0  (2.43) 
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Fig.  2.2:  Temperature  dependent  axion  potential  V{a,  T). 
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where  H(t)  =  l/2t  is  the  expansion  factor  during  the  radiation  dominated  phase 
of  the  universe.  For  rha/ma  <C  rua,  Eq.  (2.43)  has  the  approximate  solution 

(a)  ^  A{t)  cos  {via  t)  (2.44) 

where  A{t)  satisfies 

^(ma  A^)  =  -3H  {ma  A^)  .  (2.45) 

at 

Equation  (2.45)  can  be  integrated  to  yield 

{ma  A'^)f/{ma  ^2),-  =  {Rf/Ri)-^  .  (2.46) 

The  coherent  oscilJations  of  (a)  correspond  to  a  zero  momentum  Bose  condensate 
of  axions  which  is  diluted  by  the  expansion  of  the  universe.  The  energy  density 
(in  comoving  coordinates)  stored  in  those  oscillations  is  given  by 

/'«  =  5m^(a)2  +  (^)^im^?.^«?  (2.47) 

where  6,-,  the  initial  misalignment  aagle  of  the  axion  field  in  our  universe,  has 
uniform  probability  of  falling  within  the  interval  [— 7r,7r].  In  addition,  there  will 
be  dissipation  due  to  axion  decay  [a  — >  2  7]  and  [4a  — ^  2a]  which  turns  out  to  be 
negligible  in  comparision  with  Hubble  expansion.  The  contribution  to  the  present 
energy  density  has  been  calculated  to  be  [19,  20] 

fi„  A^  =  0.26  X  l0«^AJ„V(j^)"' '%.?  (2.48) 
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where  A200  is  the  QCD  scale  in  units  of  200  MeV.  In  an  inflationary  universe, 
n  =  1  and  constrains  the  RHS  of  Eq.  (2.48)  but  does  not  provide  a  definite  Hmit 

on  rua- 

Without  inflation  or  if  the  reheating  temperature  after  inflation  was  above 
the  PQ  scale,  a  different  scenario  [22,  23]  takes  place.  When  PQ  breaks,  the 
Higgs  field  a  =  vpne^'^^^'^'^  will  acquire  random  angles  in  causally  unconnected 
regions  of  spacetime.  Then  the  present  observable  universe  has  not  evolved  from 
a  domain  of  spatially  uniform  (a),  but  originated  from  many  regions  with  initially 
uncorrelated  {a).  The  Lagrangian  of  this  system  is  given  by 

L=-  df,a*d^'a  -  ^  {a*  a  -  v^q^         T  >  1  GeV  .  (2.49) 

Since  a  is  single  valued,  the  total  change  in  phase  around  any  closed  path  must  be 

A^  =  27rn  (2.50) 

where  n  is  an  integer,  called  the  winding  number.  For  n  ,1^  0,  it  is  not  possible  to 
shrink  the  loop  radius  to  zero  without  encountering  a  singularity  where  the  Higgs 
field  vanishes.  The  resulting  tubes  of  false  vacuum  are  called  strings.  The  core 
radius  is  given  by 

Sc^-yJ^ (2.51) 

V>^VpQ 
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and  the  field  far  away  from  the  core  has  the  asymptotic  form 

a  =  vpQC^^^  (2.52) 

where  6  is  the  polar  coordinate  about  the  string  axis.  The  energy  per  unit  length 
of  a  straight  global  string  is 

H^n^v%Q^-  I     1--^    27rrdr~27rn2i)2,^ln-  (2.53) 

where  the  first  term  is  the  contribution  from  the  core  and  R  is  the  cutoff  length 
which  is  typically  the  distance  between  neighbouring  strings.  Strings  with  n  ^  ±\ 
are  unstable  and  decay  into  elementary  ones.  Unlike  gauge  strings,  global  strings 
have  most  of  their  energy  outside  their  cores. 

As  the  universe  expands  kinks  will  come  within  the  horizon,  and  the  string 
network  starts  to  oscillate,  causing  the  emission  of  axions.  Furthermore,  strings 
will  intercommute  and  exchange  partners  leading  to  the  production  of  loops  which 
collapse  and  disappear.  Simulations  [24]  of  cosmic  string  evolution  suggest  that 
a  scaling  solution  will  be  quickly  reached  with  0(1)  strings  per  horizon  volume. 
This  implies  that  the  ratio  Pstring/Prad  ~  ln(/?/6),  and  strings  will  not  dominate 
the  energy  density  of  the  universe.  In  order  to  satisfy  the  scaling  solution,  axions 
have  to  be  radiated  at  a  rate 

ln„  =  -3i/«)n„  +  ^A,^^,  (2.64) 
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where 


J_  =  I j^Prad^  (2.55) 

Here  d?Pradldidu  is  the  radiation  spectrum  at  time  t,  Ua  is  the  number  density 
of  radiated  axions  and  p^^j  is  the  energy  density  of  radiation.  Equation  (2.54)  can 


"°«=i^X>$y  <=-^'" 


be  integrated  and  one  obtains 

r<     _d     p{t') 

IpQ    t 

where  tpn  is  the  time  of  PQ  breaking.  The  all  important  ingredient  of  this 
equation  is  the  radiation  spectrum.  One  group  of  authors  [22]  argues  that  axion 
radiation  from  strings  is  peaked  at  the  longest  possible  wavelength,  which  is  equal 
to  the  string  size.  Thus  uj{t)  ~  1/i  and  integration  of  Eq.  (2.56)  yields 

"a(<QCD)  ^  4Q"^a(<QCD)ln(^)  ,  ^{t)  ^  \  (2-57) 

where  ^qcd  ^^  ^^^  time  of  quark  confinement  and  the  axion  mass  turns  on.  Another 
group  [23]  cleums  that  the  spectrum  is  1/fc  with  cutoffs  at  1/t  and  1/6.  The  axion 
number  density  in  this  case  is 

"a(*QCD)  ^  4q  "^a(<QCD)  ,         ^i^)  =^  ;^  1"  (^)  (2-58) 

which  is  different  from  the  previous  estimate  by  a  logarithmic  factor  of  order  100. 
After  tqcDi  domain  walls  bounded  by  strings  will  form  which  decay  predominantly 
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into  gravitational  radiation.  Therefore  the  axion  energy  density  at  present  is 

(2.59) 


^string  -  J^^    if    ^i^)  -  ^    and 


^string  ^J^    if    ^^^^  ^  7  ^Kz)   ' 
In  summary,  cosmological  considerations  cannot  provide  a  strict  bound  on 
the  axion  mass.    If  there  was  no  inflation  or  if  inflation  ocurred  before  PQ  sym- 
metry breaking,  then  a  lower  hmit  on  rria  of  ~  (10~^  -  10~^)eV  exists,  whereas 
in  case  of  inflation  after  PQ,  no  conclusive  statement  can  be  made. 

Axions  as  a  Dark  Matter  Candidate 

Cosmologically  produced  axions  can,  for  certain  values  of  rria  and  ^,-,  make  up 
a  sizeable  fraction  of  the  present  energy  density  of  our  universe.  Axions  originating 
from  vacuum  misalignment  or  string  oscillations  are  born  cold  and  never  come  into 
thermal  equilibrium.  Thus  cosmological  axions  are  an  ideal  candidate  for  cold  dark 
matter  which  is  currently  popular  in  theories  of  galaxy  formation. 

Evidence  for  dark  matter  has  existed  for  a  long  time,  starting  with  the 
measurement  of  the  rotational  velocity  of  the  Coma  cluster.  These  measure- 
ments yielded  a  much  too  large  value  if  one  considers  only  the  masses  in  form 
of  stars.  Zwicky  [25]  apphed  the  virial  theorem  and  estimated  a  mass-to-hght 
ratio  of  M/L  ~  50/i  in  solar  units  which  is  about  20  times  larger  than  in  the 
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solar  neighbourhood.  He  concluded  that  most  of  the  mass  of  the  Coma  cluster 
is  invisible  or  dark.  More  recent  measurements  [26]  on  clusters  yield  an  average 
value  of  M/L  ~  400/i. 

Another  indication  for  dark  matter  comes  from  rotation  curves  of  spiral 
galaxies.  This  type  of  galaxy  makes  up  about  65%  of  all  observed  galaxies,  among 
them  the  Milky  Way.  Doppler  shift  measurements  of  the  rotational  velocity  of 
tracer  stars  versus  distance  from  the  center  yield  surprisingly  flat  rotation  curves. 
By  applying  Kepler's  law 

^^L(^)  ^  GmM{r)  ^2.60) 

one  expects  Vrot  to  fall  off  as  r"^'^  outside  of  the  luminous  component,  if  there  is 
no  dark  halo.  Here  m  is  the  mass  of  the  tracer  star  and  M{r)  is  the  mass  within 
radius  r  from  the  center  of  the  galaxy. 

Figure  2.3  shows  a  rotation  curve  of  NGC  3198  fitted  by  a  disk  and  a  halo  com- 
ponent. The  halo  model  used  was 

Phalo{r)  =  />/,a/o(0)(l  +  (9'')'^  (2-61) 

where  a  is  hnked  to  the  core  radius  and  7  «  2.  The  inferred  M/L  of  spiral  galaxies 
is  typically  ~  30  —  40  or  even  bigger  because  the  extent  of  the  halos  is  not  known. 
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Fig.  2.3:  Fit  to  observed  rotation  curve  with  exponential  disk  and  spherical 
halo  for  NGC  3198.  The  parameters  in  Eq.  (2.61)  are  a  -  1.3  kpc,  7  =  2.05, 
PhaloiO)  =  0.0063  M0PC-3.  (After  [27]) 
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Several  attempts  were  made  or  are  underway  to  measure  the  density  of  the 
universe  over  very  large  scales.  The  infrared  survey  [28]  done  by  the  IRAS  satelhte 
has  measured  the  density  to  lie  within  0.3  <,  fl  ^  1.2,  while  another  test  comes 
from  the  redshift-number  count  by  Loh  and  Spillar  [29]  yielding  0.7  <  9,  <  1.3. 
This  is  in  agreement  with  the  theoretically  preferred  value  f)  =  1  (M/L  ~  1600) 
as  predicted  by  inflationary  models. 

Therefore  a  great  deal  of  dark  matter  of  unknown  composition  is  likely  to 
exist.  An  important  constraint  on  its  composition  comes  from  primordial  nucle- 
osynthesis. From  the  relative  abundances  of  elements  produced  in  the  Big  Bang, 
mainly  ^H,  ^He,  '^Heand^Li,  one  can  estimate  [30]  the  density  of  baryonic  matter 
(  protons  and  neutrons  )  to  be  0.03  <  ^baryon  ^  0-12.  For  this  reason,  most 
of  the  dark  matter  has  to  be  of  some  exotic  form,  if  one  assumes  fi  =  1  as  in 
inflationary  models.  Many  candidates  have  been  put  forth,  among  them  massive 
neutrinos,  supersymmetric  particles  and  the  axion.  They  can  be  classified  into 
'cold  dark  matter'  aiid  'hot  dark  matter'.  The  former  type  {e.g.  axions)  is  slow 
moving  at  the  time  of  structure  formation  and  can  clump  on  all  scales  causing 
small  structures  to  appear  first  ('bottom-up').  The  second  group  (e.g.  hght  neu- 
trinos) is  freely  streaming  on  galactic  scales  and  produces  large  structures  first 
('top-down'). 
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It  might  be  the  case  that  neither  of  the  two  groups  can  by  itself  explain 
the  structure  of  the  universe  without  some  further  ingredient.  Axions  however 
remEiin  on  the  top  of  the  list  as  potential  candidates  for  dark  matter  because  their 
existence  is  needed  in  particle  physics. 

Detection  of  Galactic  Axions 

Turner  [20]  has  calculated  the  local  halo  density,  using  a  similar  model  as  in 
Eq.  (2.61),  to  be 

Phaio  ^  5  X  10"^^  g  cm"^  ~  0.3  GeV  cm"^  .  (2.62) 

The  halo  is  assumed  to  form  an  isothermal  sphere  with  a  Maxwellian  velocity 
distribution 

p(i;)d\  =  iVe-2^'^d\  (2.63) 

where  t3  =  3/{v  ),  {v  )^  ~  280  km  sec"  and  A'^  is  a  normalization  constant.  This 
distribution  is  broadened  by  the  earth's  rotation  around  the  center  of  our  galaxy. 
The  kinetic  energy  of  an  axion  is 

E=^ma{y-v^f  (2.64) 

where  v^  —  220  km  see"  .  The  probability  distribution  of  E  becomes 

p (E)  =  N'  e-^ ^l"^"  sinh  (/?  Ve  y/2E/ma)  (2.65) 
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which  is  plotted  in  Fig.  2.4.  The  figure  shows  that  galactic  axions  have  a  fractional 
width  l/^«~  0(10-^). 

At  present,  the  cavity  detector  is  the  most  promising  detector  for  galactic 
axions  as  proposed  by  Sikivie  [6].  It  uses  the  electromagnetic  coupling  of  the  axion, 
which  can  be  written  as 

Lay-r  =  -\  9ajj  a F^^F'^^  =  ga-n  « E  •  B  .  (2.66) 

This  coupling  allows  for  axion-photon  transitions  in  a  background  field  Bo{Eq). 
Bq  provides  a  virtual  photon  (see  Fig.  2.5)  and  the  axion  energy  is  carried  away 
by  a  single  real  photon  of  frequency 

u;  =  ma(l  +  O(10-^)).  (2.67) 

For  axion  masses  with  rua  =  O  (10~^)  eV,  the  photons  are  in  the  microwave  region, 
and  as  a  result,  the  most  suitable  detector  is  a  microwave  cavity. 

The  axion  field  can  be  described  as  a  classical  field  because  of  its  high  number 
density  (~   lO^^cm"^)  and  is  coherent  on  laboratory  scales  with  a  de  BrogUe 

wavelength  of 

,„  =  ^.10m(i5:^).  (2.68) 

ma  V  V      TUa      ' 

With  Bo  =  z  Bq,  the  coupling  of  the  axion  field  to  a  given  cavity  mode  is 

gay-rdBo  J  dhEnip{y^,t)-z  (2.69) 
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Fig.  2.4:  Kinetic  energy  distribution  of  galactic  axions  with  earth  motion 


included.  The  fractional  width  is  ~  10 
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Fig.  2.5:  Feynman  diagram  oi  a  —y  j  conversion  in  background  field  Bq. 


33 


where  E„/„(x,/)  is  the  electric  field  of  that  mode  and  V  is  the  cavity  volume. 
E„/p(x,<)  satisfies 

V  X  V  X  E„/p (x, t)-ne (x) ulip E„,p  (x,  t)  =  0  (2.70) 

where  we  have  allowed  the  dielectric  constant  of  the  material  within  the  cavity 
to  be  spatially  dependent.  From  Eq.  (2.69),  only  TM„/o  modes  couple  in  first 
approximation.  The  equation  of  motion  can  be  derived  from 

L  =  i  (e  (x)  e2  (x,  t)--B^  (x,  t))  -ga^ja  (t)  E  (x,  ^-B  (x,  t)  (2.71) 

yielding 

V  X  V  X  E  (x)  -  e  (x)  w^  E  (x)  =  u;^  gajj  a  (u)  Bq  (2.72) 

assuming  a  is  spatially  independent.  This  equation  can  best  be  analyzed  [31]  by 
expanding  the  fields  into  cavity  modes  e  (x)„/p  with  the  normahzation 

/   e  (x)  e  (x)n/p-e  (x)„///p/  d^x  =  V  6„„/  6///  6pp>  .  (2.73) 


Thus 


E  (x)  =  ^  Xnip  e  (x)n;p     and 
n,l,p 

Bo(x)  =  e(x)Bo   Y^  r]nlp^i^)nlp 
n,l,p 


(2.74) 


where 


Vnlp=^  J    i-e{x)nipd^x  .  (2.75) 
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For  a  particular  mode  one  obtains 

up' 
Klp  =  —^ o9af'ra{^)BQT]nl  (2.76) 

and  the  average  energy  in  the  cavity  mode  becomes 

U„„  =  vJ^l,/i  (2.77) 

For  a  cavity  of  quahty  Q,  Eq.  (2.77)  becomes 

2        2  2     /"  ^  (uj)  u>  du  „  , 

U„l,  =  S.„  1„„  VB„J  -j-—^^^—,^  -  .  (2. ,  8) 

The  power  from  axion  to  photon  conversions  P  —  to  V jQ  is  [6,  31] 


Pnlv  -  glj,  Cnlp  V  Bl  ^  Min  (Q,,p,  Qa)  (2.79) 


ruc 


provided  a(u;)  peaks  at  a;„/p,  and  we  have  defined  the  dimensionless  form  factor 
^nlp 


(/k  z-e(x)„/pd3xj 


For  a  cyliudrical  cavity  one  has 


(XOn)"' 


Cnlp=T-y2^l^^V^  (2.8i; 


where  xon  is  the  nth  zero  of  the  Bessel  function  Jq{x)  and  e  (x)  =  1.  It  is  preferable 
to  use  the  TMqio  mode  (Cqio  =  0.69  for  e  (x)  =  1)  because  the  higher  modes  have 
a  rapidly  decreasing  form  factor. 
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In  the  DFSZ  model,  P  from  Eq.  (2.79)  becomes 

-2S       f  ^  \  ( Bq\'^  ^       (  Pa  \  (       rua        \  (Qnlp\ 

P„/p  :.  4  X  10-     W  (— j  [—)    Cr^i,  ^5^  10-25  g/en,3J  l27r(3GHz);  VW) 

(2.82) 
where  the  quantity  Qn/p  is  the  loaded  quality  factor  of  the  cavity  given  by 

—  =  —  +  -^  .  (2.83) 

Ql     Qw     Qh 

Here  l/Qh  is  the  contribution  from  the  coupling  hole  and  l/Qw  is  the  contribution 
due  to  absorption  into  the  cavity  walls.  The  power  transferred  to  the  receiver  is 

Pr  =  ^P  (2.84) 

Qh 

where  we  have  dropped  the  mode  indices.  Since  usually  Qi  <  Qa,  one  can  look 
at  Qa/Ql  frequency  bins  simultaneously.  Because  of  the  low  signal  power  level,  it 
is  necessary  to  average  over  some  time.  The  sensitivity  of  a  microwave  receiver  is 
expressed  in  terms  of  its  system  noise  temperature  r„.  The  signal  to  noise  ratio 
s/n  for  an  integration  time  <  in  a  bandwidth  B  is 

i  =  ^^/B7=-^^/V^  (2.85) 

n        In  K^ln 

where  Tg  is  the  signal  temperature.  The  time  t  required  to  achieve  a  given  signal 
to  noise  value  s/n  is  therefore 
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The  search  rate  df/dt  at  which  frequency  bands  can  be  searched  is  proportional 
to  0-IQl)I{Q\IQ\)  which  is  maximized  for  a  Q^  value  of 

Qh      Qw      3  Qi 

The  search  rate  df/dt  for  the  optimum  coupling  configuration  (Eq.  2.87)  is  [31] 


/ 0-3 MHz N   /_^X2  /  J^x  2  /^x  4         (_Pa_Y  (^\^  (     f      \  (Q]^\ 
~\     year    J  \s/n)    VlO£/    Ut/  V/ia/o^    ^  ^n  >*    VlGHzAlO^; 


(2.88) 


CHAPTER  3 
DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 


Magnet  and  Cryogenic  System 

The  magnet  [32]  used  in  this  experiment  was  a  superconducting  (NbTi) 
solenoid  which  was  40  cm  long  and  possessed  a  bore  of  17.1  cm.  The  central  field 
was  8.6  T  at  a  current  of  88  A.  This  provided  an  average  field  of  7.5  T  in  a  volume 
of  about  8  hters.  The  magnet  was  charged  by  a  IPS-100  power  supply  [32]  and 
typically  run  in  persistant  mode  at  full  field.  Each  of  the  two  vapor  cooled  current 
leads  consisted  of  48  strands  of  10  mil  tungsten  wire. 

The  magnet,  cavity  and  cryoelectronics  were  housed  in  a  superinsulated 
dewar  [33],  which  could  hold  up  to  70  Hters  of  liquid  helium.  A  lambda  fridge 
[34]  lowered  the  bath  temperature  in  the  lower  part  of  the  dewar  to  2.2  K  (see 
Fig.  3.1).  The  fridge  (see  Fig.  3.2)  consisted  of  a  stainless  steel  tube  (1.3  cm  dia.), 
which  contained  a  heat  exchanger  and  a  flow  impedance. 

The  heat  exchanger  was  made  out  of  a  1  m  long  piece  of  1  mm  o.d.  Cu-Ni 
tubing,  while  the  impedance  was  made  out  of  a  3  cm  long  section  of  a  0.1  mm 
i.d.  Cu-Ni  tube.  Liquid  helium  entered  the  exchanger  through  a  hole  in  the 
tube  at  T  =  4  K,  was  precooled  and  exited  through  the  impedance.  A  vacuum 
pump  [35]  was  connected  to  the  steel  tube  and  caused  the  hehum  which  collected 
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(persistent   node) 


IHe    @   T=4K 


Fig.  3.1:  Sketch  (not  to  scale)  of  dewar  containing  magnet,  cavity  and  am- 


plifier. 
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Fig.  3.2:  Sketch  (not  to  scale)  of  A-fridge. 
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at  the  bottom  to  become  superfluid.  A  copper  plug  with  brush  acted  as  a  heat 
exchanger  with  the  surrounding  hehum  bath.  The  cooHng  power  of  the  fridge  was 
approximately  1  W.  Cooling  down  the  magnet  and  cavity  from  T  =  4K  to  T  =  2K 
took  about  a  day  after  which  the  temperature  stabihzed  at  the  A-transition  where 
the  heat  capacity  and  thermal  conductivity  of  He  becomes  large.  At  the  top  of 
the  bath,  the  temperature  was  T  =  4  K  at  atmospheric  pressure.  Due  to  the  poor 
thermal  conductivity  of  normal  fluid  "^He,  a  transition  region  with  a  temperature 
gradient  formed  between  the  top  and  the  position  of  the  fridge. 

The  thermal  conduction  losses  into  the  dewar  were  minimized  by  using  cryo- 
genic coaxial  cables  [36]  for  all  microwave  links  and  manganine  wires  for  the  tem- 
perature sensors.  The  helium  consumption  was  about  12  liters/day  when  running 
at  T  =  4  K  and  20  liters/day  when  using  the  A-fridge.  In  the  latter  case,  trans- 
fers were  made  every  two  days.  During  transfers,  the  bath  temperature  increased 
shghtly  and  data  acquisition  was  stopped  briefly.  Shortly  afterwards  the  bath 
stabilized  at  the  original  temperature  and  data  taking  resumed. 

The  Cavities 

Two  cylindrical  cavities  were  built  out  of  oxygen  free  copper  [37],  each  con- 
sisting of  a  tube  and  two  end  plates.    After  rough  machining,  the  surfaces  were 
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smoothed  with  #  600  sandpaper  and  finally  electropolished  [38]  to  reduce  further 
their  roughness.  The  pieces  were  then  clamped  tightly  together  in  order  to  make 
good  electrical  contact.  No  vacuum  seal  was  required  since  during  detector  oper- 
ation, the  cavity  was  filled  with  liquid  helium  which  entered  through  the  coupling 
and  tuning  rod  holes.  Figure  3.3  shows  cavity  I,  which  was  used  in  our  first  scan. 
Cavity  II  was  identical  to  cavity  I  except  that  the  radius  was  «  10  %  smaller  than 
in  cavity  I. 

Several  holes  were  made  in  the  plates  for  coupling  and  tuning  purposes. 
The  larger  one  of  the  two  coupling  ports  was  overcoupled  and  connected  to  the 
receiver.  It  consisted  of  an  inductive  loop  (shown  in  Fig.  3.4)  made  out  of  a  piece 
of  50  ^  coax  cable.  The  coupling  strength  was  adjustable  by  varying  the  insertion 
depth  of  the  loop.  The  smaller  one  consisted  of  a  weakly  coupled  probe  (shown  in 
Fig.  3.4),  which  was  used  to  measure  the  resonance  frequency  and  quality  factor 
in  transmission. 

The  cavities  were  operated  in  the  TMqio  mode,  which  yields  the  maximum 
form  factor  C.  The  electric  and  magnetic  fields  (shown  in  Fig.  3.5)  of  that  mode 
for  a  cavity  filled  with  a  medium  with  dielectric  constant  e  are  [39] 

^     ^     '  (3.1) 
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major   port 


38.1    en 


15.2    en 


Fig.  3.3:  Cavity  I  with  2  dielectric  tuning  rods.    The  larger  rod  is  moved 
laterally  while  the  smaller  one  moves  vertically. 
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Fig.  3.4:  Coupling  loop  and  coupling  probe. 
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where  Jn  is  the  Bessel  function  of  order  n,  R  is  the  cavity  radius  and 

2.405  c 
.  =  _^-  (3.2) 

is  the  resonance  frequency. 

Furthermore,  the  quahty  factor  of  that  mode  is  [39] 


Q  =  -^-  (3.3) 

^      R+L 6  ^      ' 


where  L  is  the  length  of  the  cavity  and  8  is  the  skin  depth.  At  hquid  hehum 
temperatures,  copper  is  in  the  extreme  anomalous  skin  depth  regime  such  that 
140] 

<=(#^^)"'.  (3.4) 

Here,  mg  =  9.11  x  10~^^  g  is  the  electron  mass,  e  =  4.80  x  10~^^  esu  the  electron 
charge,  vp  the  Fermi  velocity  and  n  the  conduction  electron  density.  For  Cu  [41], 
vp  =  1.57  X  10^  cm  sec~^,  n  =  8.50  x  10^^  cm~^  and  the  anomalous  skin  depth  is 

6  =  2.8  x  10~^  cm  (^— ^  j        .  (3.5) 

For  cavity  I,  the  theoretical  Q  is 

Q^^U  =  2.7  X  10^  (3.6) 
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Fig.  3.5:  Electric  and  magnetic  field  profiles  for  the  TMqio  mode. 
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while  the  measured  Q  was 

Q^all  (measured)  =  1.6  x  10^  ,  (3.7) 

which  is  in  reasonable  agreement  with  the  theoretical  value. 

Tuning  of  the  cavities  was  accomplished  by  moving  dielectric  or  metal  rods 
inside  of  the  cavity.  Dielectric  rods,  in  general,  decrease  the  resonance  frequency 
whereas  metal  rods  increase  it.  Our  dielectric  was  the  low  loss  ceramic  compound 
Mg2Ti03/MgAl204  (SMAT-9.5)  [42]  with  a  dielectric  constant  of  e  ~  9.5  and 
a  loss  tangent  tan  (5  <,  10~^.  The  loss  tangent  was  measured  by  placing  a  disk 
shaped  sample  inside  a  lead  plated  superconducting  cavity.  The  microwave  losses, 
l/Q,  are 

l/Q  =  l/Qwall  +  '^/Q dielectric  +  '^/Q coupling  ■  (3-8) 

For  a  weakly  coupled  cavity  possessing  neghgible  wall  losses,  l/Q  is  then  deter- 
mined by  the  losses  in  the  dielectric  and  tanS  -  (/Q,  where  ^s  a  geometrical 
filling  factor  of  order  unity. 

Each  cavity  contained  a  large  dielectric  or  copper  rod  which  was  moved  side- 
ways for  coarse  tuning  and  a  small  dielectric  rod  for  fine  tuning.  This  configuration 
was  chosen  in  order  to  avoid  longitudinal  mode  localization  [43],  which  leads  to 
degradation  of  the  form  factor  C. 
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The  larger  rod  was  moved  through  an  arc  shaped  slot  cut  in  both  end  caps. 
The  position  of  the  rod  was  adjusted  manually  by  rotating  a  stainless  steel  rod 
on  top  of  the  dewar.  The  movement  was  translated  by  gears  into  the  motion  of 
two  arms  supporting  the  rod  through  teflon  holders  which  were  heat  shrunk  to 
the  dielectric  (copper)  rod.  The  rod  was  further  secured  against  vibrations  by 
spring-loading  the  supports  with  BeCu  contact  strips  [44]. 

The  fine  tuning  rod  was  moved  vertically  through  a  hole  in  the  top  plate 
of  the  cavity.  The  position  of  the  rod  was  adjusted  by  a  dove  tail  slider  on  top 
of  the  cavity  and  was  driven  by  a  computer  controlled  stepper  motor  [45]  on  top 
of  the  cryostat.  The  total  rod  travel  distance  of  the  fine  tuning  rod  was  16  cm 
corresponding  to  roughly  18000  steps  of  the  motor.  The  tuning  range  of  this  rod 
was  approximately  20  MHz. 

Table  3.1  summarizes  the  parameters  of  the  cavities  used  in  this  experiment. 
Figures  3.6,  3.7  and  3.8  show  tuning  curves  calculated  for  the  case  of  a  single 
dielectric  or  metal  rod  moving  radially  inside  a  cylindrical  cavity.  The  modes  were 
obtained  by  solving  the  2D  wave  equation  using  a  relaxation  method.  Figure  3.9 
shows  the  calculated  form  factor  C.  In  calculating  resonance  frequencies  and  form 
factors,  we  neglected  the  presence  of  the  fine  tuning  rod  in  the  cavity.  Figures 
3.10  and  3.11  are  3D-plots  of  the  electric  field  of  the  TMqio  mode  for  a  cylindrical 
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Fig.  3.6:  Cavity  frequency  versus  radial  displacement  of  a  dielectric  rod  for 
R  =  7.62  cm,  r  =  0.64  cm,  L  =  38.1  cm  and  e  =  9.5.  The  TMqio  mode  used  in 
the  search  is  shown  along  with  the  TE  modes  in  its  vicinity. 
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Fig.  3.7:  Cavity  frequency  versus  radial  displacement  of  a  dielectric  rod  for 
R  =  6.83  cm,  r  =  0.64  cm,  L  =  38.1  cm  and  e  =  9.5.  The  TMqio  mode  used  in 
the  search  is  shown  along  with  the  TE  modes  in  its  vicinity. 
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Fig.  3.8:  Cavity  frequency  versus  radial  displacement  of  a  metal  rod  for 
R  =  7.62  cm,  r  =  0.62  cm  and  L  =  38.1  cm.  The  TMqio  mode  used  in  the  search 
is  shown  along  with  the  TE  modes  in  its  vicinity.  To  avoid  cluttering,  the  higher 
TM  modes  as  well  as  the  TEM  modes  are  not  drawn. 
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Fig.  3.9:  Form  factor  C  of  the  TMqiq  modes  used  in  the  search.    Curve  1 
is  calculated  for  a  radially  moving  dielectric  of  radius  r  =  0.64  cm  and  dielectric 
constant  e  =  9.5.    Cavity  radius  R  =  6.83  cm.    Curve  2  shows  the  same  for 
R  =  7.62  cm.    Curve  3  is  calculated  for  a  metal  rod  of  radius  r  =  0.62  cm  and 
cavity  radius  R  —  7.62  cm. 
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Fig.  3.10:  3D-plot  of  the  electric  field  for  the  TMqio  mode  of  a  cylindrical 
cavity  with  a  radially  displaced  dielectric  rod.  R  =  6.83  cm,  r  =  0.64  cm,  e  =  9.5 
and  rod  displacement  from  center  is  3.2  cm. 
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Fig.  3.11:  3D-plot  of  the  electric  field  for  the  TMqio  mode  of  a  cyhndrical 
cavity  with  a  radially  displaced  metal  rod.  R  -  7.62  cm,  r  =  0.62  cm  and  rod 
displacement  from  center  is  2.1  cm. 
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Fig.  3.12:  Measured  mode  crossing  between  the  TE112  mode  and  the  TMqio 


mode. 
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Fig.  3.13:  Transmitted  power  through  cavity  II  near  the  TMqio  resonance. 
The  cavity  was  overcoupled  with  Ql  =  Qwalll"^  ~  40000. 
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Fig.  3.14:  Reflected  power  off"  major  port  of  cavity  II  near  the  TMqio  reso- 
nance. The  return  loss  is  r;  10.5  dB. 
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Table  3.1    Summary  of  cavity  parameters.  L  =  38.1  cm  in  all  cases. 


R  [cm] 

tuning  rod 

tuning  range  [GHz] 

Qwall 

C 

7.62 
6.83 
7.62 

dielectric 

dielectric 

metal 

1.32-1.44 
1.44-1.60 
1.56-1.85 

1.5  X  10^ 
1.2  X  10^ 
0.5  X  10^ 

0.50 
0.45 
0.60 

cavity  with  a  displaced  dielectric  and  metal  rod  respectively.  As  seen  from  the 
tuning  curves  in  Figs.  3.6,  3.7  and  3.8,  a  nvimber  of  crossings  between  the  TMqio 
mode  and  TE  modes  occur.  In  the  neighbourhood  of  the  crossings,  the  modes  will 
mix  and  repel  each  other  leading  to  holes  in  the  spectrum.  Figure  3.12  displays 
such  a  crossing  as  measured  in  cavity  I.  The  holes  are  typically  a  few  hundred 
kHz  wide.  Figure  3.13  shows  a  transmission  spectrum  of  the  overcoupled  cavity  H 
near  the  TMqio  resonance,  with  Q^all  =  ^Qhole^  while  Fig.  3.14  shows  the  power 
reflected  off  the  major  port  for  the  same  resonance. 


The  Microwave  Amplifiers 

The  cryogenic  microwave  amplifiers  were  purchased  from  Berkshire  Tech- 
nologies [46].  The  first  of  the  three  stages  consisted  of  a  high-electron-mobility 
transistor  (HEMT),  while  the  remaining  stages  utilized  ordinary  GaAs  microwave 
transistors.  Upon  cooling,  the  noise  temperature  of  these  devices  drops  dramati- 
cally. A  photograph  of  one  of  the  amplifiers  is  shown  in  Fig.  3.15,  whereas  Fig.  3.16 
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shows  an  amplifier  located  in  the  cryostat. 

The  noise  temperature  of  the  amplifiers  was  carefully  measured  using  a  vari- 
able temperature  technique  [47].  A  50  0  termination  acted  as  the  noise  source 
and  was  connected  to  the  amplifier  input.  The  amplifier  excess  noise  is  obtained 
by  varying  the  temperature  of  the  load  and  measuring  the  output  noise  of  the 
amplifier.  Figure  3.17  shows  the  test  setup  including  load  and  amplifier. 

The  termination  consisted  of  a  metal  film  rod  resistor  [48]  on  a  BeO  substrate 
with  one  end  soldered  to  the  center  conductor  of  a  section  of  cryogenic  coax  cable. 
The  other  end  was  soldered  to  a  copper  plate,  which  had  a  thermometer  and  heater 
attached  to  it.  The  time  constant  of  a  few  seconds  was  determined  by  a  thermal 
link  consisting  of  two  10  cm  long  strands  of  AWG-18  copper  wire  between  the 
copper  plate  and  the  helium  bath.  The  microwave  feedthrough  of  the  vacuum  can 
was  a  commercial  heiinetir  connector  [49]  with  the  rubber  o-ring  replaced  by  an 
indium  o-ring.  The  emitted  noise  power  of  the  termination  in  a  bandwidth  B  is 

P  =  kBTB{\-\p\^)  (3.9) 

where  k^  is  Boltzmann's  constant,  T  is  the  load  temperature  and  p  is  the  reflection 
coefficient.  The  power  return  loss  of  the  unit  was  measured  to  be  greater  than 
20  dB,  thus  giving  only  a  small  correction  to  the  perfectly  matched  case.    The 
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Fig.  3.15:  Photograph  of  L-band  HEMT  amplifier  with  top  lid  removed.  The 
input  is  to  the  left  and  three  amplification  stages  are  visible. 
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Fig.  3.16:  Photograph  of  cryostat  showing  magnet  and  cavity.  The  fine 
tuning  drive  is  in  the  center,  and  the  adjustable  couphng  is  to  the  right.  The 
ampHfier  and  circulator  are  located  near  the  top. 
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Fig.  3.17:  Setup  of  noise  temperature  measurement. 
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electrical  loss  of  the  line  between  load  and  amplifier  was  less  than  0.1  dB  and  was 
neglected. 

The  noise  power  at  the  output  of  the  amplifier  is 

P  =  kB{Ta  +  T)BG  (3.10) 

where  Ta  nad  G  are  the  amplifier  noise  temperature  and  gain  repectively.  This 
power  was  further  amplified  by  a  room  temperature  amplifier  [50]  and  finally  de- 
tected by  a  HP  8569B  spectrum  analyzer  [51].  For  a  given  frequency,  the  noise 
temperature  is  obtained  by  measuring  P  (Eq.  3.10)  for  different  T's  and  fitting 
these  data  with  a  straight  line  to  extract  Tq.  The  results  for  two  cryogenic  ampli- 
fiers are  given  in  Figs.  3.18  and  3.19. 

For  these  ultra  low  noise  amplifiers  it  is  necessary  to  take  into  account  the 
noise  contributions  from  the  other  stages  in  the  chain.  The  formula  for  cascaded 
amplifiers  is 

T  =  Ti  +  T2/G1  +  T3/G2  +  •••  (3.11) 

where  the  Tn's  and  Gn's  are  the  noise  temperatures  and  gains  respectively.  The 
post  amplifier  had  a  gain  of  50  dB  and  a  noise  temperature  of  200  K,  making 
all  following  contributions  negligible.  The  two  cryogenic  amplifiers  had  gains  of 
34  dB  and  35  dB  with  minimum  noise  temperatures  of  order  Tq  «  3  K.  Thus  the 
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Fig.  3.18:  Noise  temperature  of  Berkshire  L-1.5-30H  at  a  physical  tempera- 
ture of  r=  4  K.  The  power  gain  is  ~  34  dB  over  the  shown  band. 
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Fig.  3.19:  Noise  temperature  of  Berkshire  S-2.3-30HR  at  a  physical  temper- 
ature of  T=  4  K.  The  power  gain  is  ~  35  dB  over  the  shown  band. 
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noise  from  the  post  amplifier  adds  a  small  amount  of  extra  noise  to  the  system 
( Ri  0.1  K)  which  is  included  in  Figs.  3.18  and  3.19. 

In  the  experiment,  magnetically  shielded  cryogenic  circulators  [52]  were  in- 
serted between  cavity  and  amplifier  for  impedance  matching.  This  procedure  was 
necessary  because  the  cavity  impedance  changed  rapidly  across  the  resonance,  and 
the  amplifier  was  optimized  for  a  constant  input  impedance  of  50  ft.  The  circu- 
lator uses  its  third  arm  to  absorb  any  reflected  power  coming  from  the  amplifier 
input.  The  third  port  of  the  circulator  was  also  used  to  bring  microwave  power  to 
the  major  port  for  measurement  of  its  reflection  coefficient.  In  order  to  keep  room 
temperature  noise  from  entering  the  amplifier  or  cavity,  a  50  dB  attenuator  [53] 
was  added.  Assuming  the  coax  cable  is  properly  terminated  at  room  temperature, 
the  noise  power  per  unit  bandwidth  emitted  by  the  attenuator  is  [54] 

kBT  =  kBToa  +  kBTattil-a)  "  (3.12) 

where  a  is  the  power  loss  and  Tq  ~  300  K.  Therefore,  the  attenuation  is  sufficient 
to  reduce  the  incoming  noise  to  a  level  equivalent  to  the  bath  temperature. 
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Room  Temperature  Data  Acquisition  Hardware 

Figure  3.20  shows  the  complete  experimental  setup  and  data  acquisition 
system.  The  output  of  the  cryogenic  amplifier  was  fed  into  a  post-amplifier  [50] 
and  subsequently  into  the  first  mixing  stage. 

The  image  reject  mixer  (see  Fig.  3.21)  was  assembled  from  commercially  available 
components.  The  principle  of  an  image  reject  mixer  is  to  suppress  the  unwanted 
sideband  through  destructive  interference.  The  image  reject  ratio  of  the  mixer 
refers  to  the  power  ratio  at  the  mixer  output  of  the  two  sideband  frequencies. 

For  a  given  intermediate  frequency  wjp,  the  mixer  is  sensitive  to  two  side- 
bands represented  by  the  voltages  V  =  cos[(u;  ±wip)<],  where  u  is  the  LO  fre- 
quency. The  RF  hybrid  [55]  splits  this  signal  into  two  components  with  a  90° 
phase  shift  relative  to  each  other: 

cos    (u;  ±a;ip)f    ->  — p  I         .,      ,         .^,       /oil-  (o.lcJj 

It  is  followed  by  two  identical  doubly  balanced  mixers  |56]  which  down  converted 

the  signal  to  an  IF  frequency,  i.e. 

_]_^cos[{u±cJiY)t]  \        J_  /cos[wiF<]  \  .3J4^ 

y/2\cos[{uJ±u;iY)t  +  7T/2]J        y/2  \cos[uJiYt  ±  n/2]J  ' 

where  the  mixer  losses  have  been  neglected.  The  LO  [57]  output  power  level  was 

set  to  10  dBm  and  was  divided  in  phase  by  a  power  splitter  [58].  The  IF  frequency 
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Fig.  3.20:  Block  diagram  of  experimental  setup.  The  room  temperature  part 
consists  of  a  two-stage  superheterodyne  receiver  followed  by  a  real  time  spectrum 
analyzer. 
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Fig.  3.21:  Block  diagram  of  image  reject  mixer.  The  image  reject  ratio  was 
Z  20  dB  in  the  1-2  GHz  band. 
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of  10.7  MHz  was  determined  by  the  center  frequency  of  the  crystal  filter  following 

the  mixer.   Finally  the  IF  signals  were  combined  in  an  IF  hybrid  [59]  where  the 

the  two  sidebands  appeared  at  two  different  ports.    The  effect  of  the  hybrid  can 

be  represented  by 

+    /cos[u;ipt]' 


J_  fcos[uiYt]  \         V         0         /  ,^^^. 

V^^cos[a;iFt±7r/2]y'    -    /  0  \     ' 

\cos[coiYt  +  TT/2]J 

where  the  upper  and  lower  entries  in  the  bracketed  expressions  on  the  RHS  give 
the  voltages  at  the  two  output  ports.  As  a  result,  the  two  sidebands  appear  at 
different  ports  and  the  undesired  one  is  terminated.  The  RF  loss  of  the  unit  was 
PS  8  dB  and  the  image  reject  ratio  was  measured  to  be  >  20  dB  (see  Fig.  3.22). 

The  signal  was  further  amplified  by  two  IF  amplifiers  [60],  with  a  combined 
gain  of  60  dB.  The  signal  then  passed  through  a  8  pole  crystal  filter  [61]  before 
entering  the  second  mixing  stage.  The  filter  had  a  center  frequency  of  10.7  MHz 
and  a  3  dB  pass  band  of  30  kHz.  The  attenuation  was  w  3  dB  over  the  useful  band 
(see  Fig.  3.23)  and  was  sHghtly  temperature  dependent.  To  stabiHze  the  response, 
the  filter  was  enclosed  in  a  box  and  its  temperature  held  at  T  =  (29  ±  .1)  °C  by  a 
temperature  controller  [62]. 

The  second  mixing  stage  comprised  a  single  double  balanced  mixer  [63],  and 
a  local  oscillator  [64].  The  audio  signal  covering  DC  to  ss  30  kHz  was  finally  passed 
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Fig.  3.22:  Image  reject  ratio  versus  frequency  of  first  stage  mixer  for  an  IF 
of  10.7  MHz.  The  solid  curve  is  for  /lo  >  /signah  while  the  dotted  curve  is  for  the 
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Fig.  3.23:  Response  of  crystal  filter. 
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to  an  AC  amplifier  and  sampled  by  an  A/D  converter  with  a  resolution  of  16  bit. 
The  samphng  rate  was  set  to  70  kHz  in  order  to  oversample  slightly.  The  signal 
was  then  spectrum  analyzed  and  averaged  in  real  time.  For  each  spectrum,  64 
samples  were  taken  yielding  a  32  point  power  spectrum.  The  samples  were  stored 
in  a  buffer  memory,  and  while  one  set  of  samples  was  taken,  the  previous  array 
was  fourier  analyzed  and  averaged.  The  power  spectra  were  computed  by  a  TI 
TMS320C25  digital  signal  processing  chip  which  is  a  fixed  point  device  with  an 
instruction  cycle  of  100  ns.  All  digital  processing  components  including  a  RAM 
of  64  kbytes  for  data  and  program  storage  were  on  a  single  board  [65]  which  was 
integrated  into  a  Z-386  workstation  [64]. 


CHAPTER  4 
DATA  ACQUISITION  AND  ANALYSIS 

The  Search  Process 

Data  taking  in  this  experiment  was  a  time  consuming  and  repetitive  task. 
Therefore  as  much  as  possible  of  the  process  was  automated  and  runs  were  left 
unattended  for  many  hours. 

The  main  data  acquisition  program  was  called  FNDAX2  and  is  listed  in 

the  Appendix.  It  controls  various  intruments  with  GPIB  commands  or  by  RS232 

serial  line.  Before  the  program  is  run,  the  coupling  strength  of  the  major  port  is 

set  manually  to  the  desired  value.  The  program  starts  by  measuring  the  resonance 

frequency,  Q  and  reflection  coefficent  of  the  TMqiq  mode.  For  this  purpose,  a  set 

of  three  HP33311B  microwave  switches  [51]  were  installed  and  controlled  by  a 

HP59306A  relay  actuator  [51].    First  CW  power  from  the  LO  is  directed  to  the 

weakly  coupled  minor  port  and  the  transmitted  power  is  measured  for  several 

frequencies  around  the  resonance  by  a  HP438A  power  meter  [51].  These  data  are 

then  fitted  by  a  Lorentzian  to  obtain  resonance  frequency  and  Q.  The  LO  is  then 

switched  to  the  cryogenic  circulator  and  reflected  power  is  measured  on  and  off 

resonance  which  yields  the  coupling  strength  of  the  major  port. 

At  this  stage,  the  LO  enters  the  data  acquisition  mode  and  is  switched  to  the 
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image  reject  mixer.  The  host  PC  then  starts  the  spectrum  analysis  program  which 
was  downloaded  into  the  DSP  board  prior  to  the  run.  The  A/D  converter  samples 
in  interrupt  mode  the  audio  signal  coming  from  the  amplifier  chain  and  stores 
the  data  in  one  of  two  buffer  memories.  The  TMS320C25  chip  in  the  meantime 
works  on  the  other  buffer  and  computes  a  64  complex  point  FFT  resulting  in  a 
32  point  power  spectrum.  The  time  for  one  spectrum  is  0.9  ms  and  a  total  of 
10  spectra  are  taken  and  added  together  in  approximately  90  seconds.  After  the 
spectrum  is  completed,  the  host  uploads  the  data  and  normalizes  it.  In  order  to 
get  a  fiat  baseline,  the  frequency  response  of  the  crystal  filter  has  to  be  divided 
out.  Moreover,  the  noise  output  of  the  cryogenic  HEMT  amplifier  can  deviate 
from  a  flat  spectrum  if  temperature  imbalances  in  the  circulator-amplifier-cavity 
system  exist.  To  account  for  all  those  effects,  a  running  average  is  calculated  from 
the  last  30-40  scans  and  used  for  normalization. 

For  each  bin,  the  deviation  from  the  mean  is  calculated  according  to 

#  Sigma  =:^^V7V  (4.1) 

where  N  is  the  number  of  averages  and  the  spectrum  is  displayed  on  the  screen. 
The  program  searches  the  spectrum  for  peaJcs  above  2  a  in  single  bins  and  com- 
binations of  two  neighbouring  bins.  The  bin  width  of  1  kHz  roughly  matches  the 
expected  width  of  the  axion  line.  If  a  peak  appeared,  another  spectrum  is  taken 
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and  compared  with  the  first.  In  the  first  cavity  scans,  the  spectra  were  added  to- 
gether, but  in  the  later  scans  the  two  spectra  were  checked  for  coincidence  peaks. 
This  process  was  repeated  up  to  four  times  and  the  peak  was  flagged  for  later 
rechecks  if  it  proved  to  be  statistically  significant.  Each  spectrum  is  written  into 
a  file  along  with  the  values  of  the  cavity  parameters  and  the  temperature  of  the 
bath.  The  computer  then  resets  the  cavity  frequency  by  sending  a  number  of 
pulses  to  the  stepping  motor  and  the  cycle  repeats  itself.  A  typical  run  produces 
about  400  spectra  corresponding  to  a  tuning  range  of  about  8  MHz.  The  data  are 
later  backed  up  on  floppy  disks  for  permanent  storage. 

Occasionally,  a  test  signal  was  introduced  into  the  system.  The  CW  signal 
was  generated  by  a  HP8350A  [51]  sweeper  phase  locked  to  an  EIP578  [67]  mi- 
crowave counter  for  frequency  stability.  The  power  was  attenuated  to  ~  10~  W 
and  fed  into  the  minor  port  of  the  cavity.  This  procedure  served  mainly  for  the 
frequency  calibration  of  the  detector. 

The  Data 

To  this  date  (July  90),  five  scans  have  been  completed.  Cavity  I  with  a 
dielectric  tuning  rod  was  run  twice  at  T  =  2  K  and  T  =  4  K,  while  cavity  II  and 
cavity  I  with  a  metal  rod  were  run  once  at  T  =  2  K. 
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A  typical  spectrum  obtained  during  a  scan  is  shown  in  Fig.  4.1  whereas 
Fig.  4.2  is  an  example  of  a  spectrum  with  a  candidate  pealc.  A  nimiber  of  similar 
peaks  were  found  and  rechecked  after  the  scan  with  the  magnetic  field  turned  off. 
None  of  the  candidate  peaks  survived  this  test,  i.e.,  all  persisted  when  the  field 
was  off.  Most  of  the  peaks  were  narrow  (<  200  Hz)  and  are  assumed  to  be  pickup 
from  computer  clocks  in  the  vicinity  of  the  experiment.  The  number  of  peaks  were 
quite  high  and  could  have  been  reduced  by  putting  the  apparatus  into  a  screened 
room.  The  scans  have  some  holes  in  the  covered  frequency  range  due  to  crossings 
with  TE  modes.  The  data  from  the  four  scans  is  summarized  in  Table  4.1. 

Table  4.1 


Scan 

Frequency  Coverage  [GHz] 

Holes  [GHz] 

Peaks 

I,II 

1.32-1.44 

1.346096-1.346382 

32 

III 

1.44-1.60 

1.444693-1.445210 
1.449171-1.450566 

62 

IV 

1.60-1.63 

1.602404-1.602991 

18 

V 

1.80-1.83 

- 

2 
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Frequency  (kHz)  -  1.603021  GHz 


Fig.  4.1:  Sample  spectrum  of  cavity  output,     la  corresponds  to  a  power 
P  ~  3  X  10~^^  W  in  the  cavity. 
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Frequency  (kHz)  -  1.533251  GHz 


Fig.  4.2:  Spectrum  with  candidate  peak. 
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Limit  on  the  Electromagnetic  Coupling  of  the  Axion 

Since  no  peak  with  the  signature  of  an  axion  signal  was  found,  it  is  possible 
to  put  an  upper  limit  on  the  electromagnetic  coupling  of  the  axion,  provided  the 
galactic  halo  is  axionic.  The  noise  fluctuations  in  the  detector  have  a  Gaussian  dis- 
tribution, and  one  can  calculate  the  probability  for  a  given  signal  to  be  detectable. 
Our  detector  is  looking  for  2<7  peaks  above  the  noise  floor.  The  probability  for 
a  negative  fluctuation  of  >  2a  is  0.025  or  in  other  words,  a.  >  Aa  signal  will  be 
seen  at  a  97.5%  confidence  level.  In  terms  of  the  system  noise  temperature  Tn, 
bin  bandwidth  B  and  #  of  averages  N,  the  minimum  detectable  signal  power  is 
(at  97.5%  C.L.) 

P^i^  =  4kBTnB/VN  (4.2) 

where  kg  is  Boltzmann's  constant.  Tn  is  given  by 

Tn  =  Tf,ath  +  Tamp  (4.3) 

where  Tamp  is  almost  solely  determined  by  the  contribution  from  the  cryogenic 
amplifier.  Table  4.2  summarizes  the  parameters  entering  the  calculation  of  Pmin 
for  the  different  runs. 

Scans  I  and  II  are  special  in  that  the  same  frequency  range  was  covered 
twice.  By  calculating  the  combined  probability  for  finding  a  peak,  one  obtains  a 
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Table  4.2 


Scan 

N 

B  [kHz] 

Tbath  [K] 

Pmin  [W] 

I,II 

lo^lo^ 

1.1 

4.2,2.2 

0.7  X  10-21 

III 

10^ 

1.1 

2.2 

1.0  X  10-21 

IV 

2x  10^ 

1.1 

2.2 

1.3  X  10-21 

V 

10^ 

1.1 

2.2 

1.1  X  10-21 

power  somewhat  smaller  than  in  Eq.  4.2.  The  values  for  Pmin  ^^  Table  4.2  are 
to  be  compared  with  the  power  from  axion  to  photon  conversion  going  into  the 
detector.  With  the  cavity  overcoupled  {Qu>all  —  ^Qhole)i  oiily  ^  fraction  2/3  of 
the  power  in  Eq.  (2.82)  is  leaving  the  cavity.  Moreover,  the  cavity  response  to 
axion  decay  is  falling  off  away  from  resonance,  and  the  power  is  further  reduced 
by  a  factor  a  with  0.5  <  a  <  1.0. 

Table  4.3  lists  the  average  power  from  a  — +  7  conversion  for  the  individual 
scans.  Also  shown  is  the  axion  mass  and  the  upper  limit  on  g^^j  (for  pa  =  Phalo)^ 
which  is  plotted  in  Fig-  4.3.  This  limit  is  valid  assuming  the  power  P^  is  falling 
in  a  single  bin  and  is  worse  by  a  factor  Vz   if  it  falls  into  two  bins. 


Table  4.3 
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Scan 

Pd  [w] 

ma  [10-6  eV] 

gl^^  [GeV-2] 

1,11 

1.3  X  10-24 

5.46-5.95 

<  3.3  X  10-28 

III 

8.3  X  10-25 

5.95-6.62 

<  9.0  X  10-28 

IV 

8.4  X  10-25 

6.62-6.77 

<  1.3  X  10-2^ 

V 

1.0  X  10-24 

7.47-7.60 

<  1.2  X  10-27 
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Fig.  4.3:  Experimental  limit  on  the  electromagnetic  coupling  gajj  of  the 
axion.  Also  shown  is  the  limit  obtained  by  the  RBF  collaboration  [68].  The  solid 
and  dashed  straight  lines  are  the  theoretical  values  of  ^077  for  the  DFSZ  axion 
and  the  KSVZ  axion  (E/N  =  0). 


CHAPTER  5 
SUMMARY  AND  CONCLUSION 


The  experiment  described  herein  represents  one  of  the  first  attempts  to  de- 
tect dark  matter  in  our  galaxy.  The  precise  nature  of  the  dark  matter  is  not 
known,  but  axions  are  a  serious  candidate  for  it.  The  first  cosmic  axion  search 
was  carried  out  in  a  pioneering  experiment  by  the  RBF  collaboration  [68].  This 
experiment  is  the  second  of  its  kind. 

In  the  frequency  range  covered,  no  axion  signal  was  found.    Our  detector 

was  probably  not  sensitive  enough  to  detect  axions,  even  if  they  existed  in  our 

galactic  halo  and  with  a  mass  within  the  searched  range.  To  have  a  good  chance 

of  detecting  galactic  axions,  the  sensitivity  of  the  apparatus  needs  to  be  improved 

by  roughly  three  orders  of  magnitude  in  case  of  the  DFSZ  axion  (see  Fig.   4.3). 

Nonetheless,  the  present  data  constrain  the  coupling  gayj  of  any  type  of  axion  if 

they  are  the  dark  matter  in  the  halo.  The  limits  shown  in  Fig.  4.3  are  expressed 

in  terms  of  the  experimental  sensitivity,  which  goes  like  gayj-  The  limit  on  g^yy 

obtained  by  the  RBF  collaboration  was  improved  by  roughly  an  order  of  magnitude 

over  the  frequency  band  covered  so  far.  This  improvement  was  mainly  a  result  of 

a  lower  system  noise  temperature  resulting  from  the  use  of  HEMT  amplifiers  and 

of  lower  bath  temperatures. 
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Besides  setting  a  limit  on  the  axion  coupling,  these  two  prototype  exper- 
iments have  also  helped  to  gain  know-how  for  the  construction  of  an  improved 
detector  system.  Our  microwave  receivers  have  a  demonstrated  sensitivity  to 
signals  as  small  as  lO"'^^  W.  Since  there  is  not  much  hope  to  greatly  improve 
the  amplifier  technology  in  the  immediate  future  and  since  steady-state  magnetic 
fields  are  at  present  limited  to  about  20  T,  the  only  possible  way  to  reach  the 
DFSZ  limit  in  the  near  future  is  by  increasing  the  cavity  volume. 

If  the  present  experimental  parameters  Tn,  C,  Q,  Bq  and  d(ln/)/dt  are  not 
appreciably  changed,  achieving  the  required  sensitivity  means  increasing  the  vol- 
ume to  several  thousand  hters.  The  fundamental  resonant  frequency  of  such  a 
cavity  is  around  100  MHz.  Searching  for  the  axion  at  those  frequencies  would 
be  as  interesting  as  at  1  GHz  or  10  GHz  given  the  uncertainties  in  cosmological 
calculations  of  the  axion  density  as  a  function  of  the  axion  mass  ma- 

A  realistic  search  should  of  course  cover  a  very  large  frequency  band  while 
maintaining  high  sensitivity  and  large  volume.  The  diflSculties  associated  with 
realizing  these  constraints  increase  rapidly  with  increasing  frequency  [41].  The 
most  promising  approach  at  present  seems  to  be  dividing  the  avcdlable  volume  into 
smaller  and  smaller  cells  with  growing  resonance  frequencies.  The  cavities  have 
to  be  ganged  together  to  the  same  frequency  and  their  output  combined  in  phase. 
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The  complexity  of  such  a  system  will  limit  the  number  of  cells  to  a  few  himdred 
and  allow  the  detector  to  work  up  to  a  few  GHz.  The  higher  frequency  (mass) 
ranges  could  then  be  explored  by  a  detector  with  a  much  higher  magnetic  field 
over  a  smaller  volume.  Finally,  the  chances  for  a  sucessful  axion  experiment  could 
be  greatly  improved,  if  the  allowed  window  of  axion  masses  could  be  narrowed. 


APPENDIX  A 
CAVITY  DEVELOPMENT 


One  of  the  challenges  in  the  design  of  a  full  scale  axion  cavity  detector  is 
the  development  of  large  volume  /  high  frequency  microwave  cavities  [43].  The 
ideal  cavity  should  be  operated  in  the  lowest  TM  mode  for  maximizing  the  form 
factor  C.  One  way  of  raising  the  fundamental  frequency  is  to  insert  metal  posts 
into  the  cavity.  This  would  mostly  affect  the  TM  modes,  because  the  boundary 
conditions  for  these  modes  require  £J^°"  to  vanish  on  all  metal  surfaces.  On  the 
other  hand,  the  posts  would  not,  in  a  first  approximation,  change  the  TE  modes, 
whose  boundary  conditions  require  the  normal  derivative  of  H^^  to  be  zero  on  the 
surface.  Since  the  mode  density  of  the  TE  modes  increases  as  /  ,  the  frequency 
lifted  TM  mode  would  have  to  cross  more  and  more  TE  modes  during  tuning.  This 
mode  crossing  problem  will  eventually  become  so  serious  that  the  holes  arising 
from  the  crossings  woiild  consume  a  large  portion  of  the  searched  band.  Another 
danger  associated  with  using  many  posts  is  transverse  mode  localization.  This 
eflPect  is  caused  by  asyraraetries  in  the  post  configuration  leading  to  unpredictable 
localizations  of  the  fields. 

One  alternative  to  using  posts  is  to  employ  identical  cavities  and  combine 
their  output  signals.  The  previously  described  problems  of  localization  and  reso- 
nance crowding  would  disappear  but  at  the  cost  of  larger  complexity.  The  cavities 
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would  have  to  be  tuned  individually  to  the  same  frequency  and  their  power  com- 
bined in  phase  and  brought  to  the  front  end  of  a  single  amplifier. 

We  tested  this  method  at  room  temperature  for  the  simplest  case  of  two 
cavities.  The  test  setup  is  shown  in  Fig.  A.l.  The  empty  cavities  had  TMqio 
resonance  frequencies  of  1.67  GHz  with  Quiall  ~  3  x  10  .  The  frequencies  differed 
initially  from  each  other  by  0.7%,  but  could  be  made  to  match  by  inserting  a  small 
teflon  tuning  rod  (0.4  cm  in  diameter)  into  one  of  the  cavities. 

Two  identical  Wilkinson  powerdividers/combiners  [58]  were  used  to  split  and 
combine  the  signals  in  phase.  A  schematic  circuit  diagram  is  displayed  in  Fig.  A. 2. 
A  signal  of  voltage  V  ~  e~*'^'  entering  on  the  left  of  Fig.  A. 2  is  split  into  equal 
amplitude  signals  V  ~  e~"^'/-\/2  on  the  right.  When  used  as  a  combiner  with 
amphtudes  ae~'^'  and  6e-*"'-'<^,  the  output  voltage  is  (a  +  be'"'^)  e''"^^ /y/2, 
and  power  is  dissipated  in  the  internal  resistors  unless  the  two  inputs  are  exactly 
balanced.  In  the  experiment,  one  of  the  ports  of  each  cavity  was  critically  coupled, 
while  the  other  was  weakly  coupled. 

The  cavity  can  be  modeled  by  a  series  iJLC-circuit  with  resonance  frequency 
uJq  =  \/LC  and  Qw  =  uqL/R.  The  input  impedance  of  a  single  port  cavity  near 
resonance  is  given  by 
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Fig.  A.l:  Test  setup  of  two  frequency-matched  cavities. 
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Fig.  A. 2:  Circuit  diagram  of  Wilkinson  powerdivider/combiner. 
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where  ^  =  Qw/Qh  '^^  the  coupling  parameter  and  Zq  is  the  characteristic  hne 

impedance,  usually  50  i7.  The  input  impedance  is  equal  to  the  line  impedance  for 

a  criticcdly  coupled  cavity  (^  =  1)  on  resonance,  making  the  cavity  reflectionless. 

For  a  two  port  cavity,  Eq.    (A.l)  still  applies  but  R  is  modified  to  include  the 

additional  losses  from  the  second  port.  A  signal  transmitted  through  a  critically 

coupled  cavity  (while  the  second  port  is  weakly  coupled)  will  suffer  a  phase  shift 

V  ~  e'"  with 

tana^  (-—- u;L) /2R  c:^  2Ql'^  ~  ^^  =x  (A.2) 

Va;C  /  WQ 

and  the  amplitude  is  decreased  by  a  factor  1/vl  +  a^  relative  to  the  value  on 


resonance. 


The  output  voltage  of  the  power  combiner  for  the  LO  frequency  at  the  reso- 
nance frequency  of  the  first  cavity  and  variable  resonance  frequency  of  the  second 
cavity  is  then 

and  the  power 

\V  \2      1  /  3     \ 

TT     =7(1  +  :^ 9)  •  (A-4) 

\Vq\         4\        1  +  x^J  ^       ^ 

This  relationship  is  plotted  in  Fig.  A.3.  Also  shown  are  the  measured  values 
which  are  in  good  agreement  with  the  theoretical  curve.  Operation  of  a  multi- 
cavity  detector  would  require  keeping  the  cavity  resonance  frequencies  within  a 


91 


Fig.  A. 3:  Transmitted  power  for  the  configuration  shown  in  Fig.  A.l.  The 
LO  frequency  is  fixed  at  the  resonance  frequency  of  one  cavity  while  the  second 
cavity  is  tuned  with  a  teflon  rod. 
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few  percent  of  the  cavity  width  in  order  to  avoid  mismatch  losses.  In  summary, 
this  experiment  has  demonstrated  the  feasibihtiy  of  using  phased  multicell  cavities 
as  axion  cavity  detectors. 


APPENDIX  B 
SOFTWARE  LISTINGS 


program  fndax2 


data  acquisition  program  for  axion  search 


character  f rec*20 , f ilnm*20 , f ilcal*20 , card*2 ,bell*l , dot*l , ampl*10 

character  dir*l 

integer*4  spein(32) ,speina(32) ,nav,np,npar ,ncal 

parameter  (np=8 ,ncal=40) 

external  funct 

real*8  spel(32) ,spe2(31) ,devl(32) , dev2(31) ,oraega,devs(5 , 32) 

real*8  cal(32) ,f ,df ,aa(np) ,bb(np) ,parmt(3) , del, fold, dw 

real*8  chisq,chipr  ,alamda,parmr(3)  ,pl  ,pr  ,pin,  f  1 ,  fr  ,aw,avc,pref 

real*8  stdev,call(0:ncal-l, 32) ,skk(32) , dell , f start , f stop 

logical  peakl(32) ,peak2(31) ,peaklp(32) ,peak2p(31) ,peak,peakh 

data  frec(l:2)  /'fr'/.  frec(13:14)  /'gz'/,  fsam  /70000./, 

1  stdev  /l.Od-7/,  nav  /lOOOOO/,  bell(l:l)  /7/,  tc  /4.3/, 

2  del  /2.0d-5/,  npar  /3/,  dot  /'  .' /  ,rma.y.   /4/, 

3  ampl(l:2)  /'ap'/,  ampl(8:9)  /'db'/  ,amp  /-lO.O/ 

instrument  initialization  : 
open  (1 , file-' ieeeout' , status='old' ,access—' sequential' ) 
open  (2 , f ile=' ieeein' , status='old' ,access=' sequential' , f orm= 


'binary') 

write  (1,*)  'reset' 

write  (1,*)  'remote 

01,02,11,13' 

write  (1,*)  'clear  01,02,11,13' 

write  (1,*)  'output 

01 ; %A1B2A3 ' 

write  (1,*)  'output 

02 ; amOpc ' 

write  (1,*)  'output 

11;  *  f4  r4' 

write  (1,*)  'output 

13;kb96en' 

rewind  1 

write  (*,10) 

read  (*,*)  f 

write  (*,12) 

read  (*,*)  parmt(2) 

write  (*,9)  amp 

read  (*,*)  amp 

write  (ampl(3:7) ,8) 

amp 

write  (*,713) 

read  (*,'(a)')  dir 

if  (dir.eq. 'I 

.or. dir.eq. ' 

then 
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mult-  -1 
else 

mult-  1 
endif 
write  (*,13) 
read  (*,*)  nstep 
write  (*,14) 
read  (*,'(a)')  filnm 
ilen-  0 
15    ilen-  ilen  +  1 

if  (f ilnm(ilen+l: ilen+1) .ne . '  ')  goto  15 
write  (*,16) 
read  (*,*)  ns 
write  (*,716)  nmax 
read  (*,*)  nmax 
write  (*,17)  tc 
read  (*,*)  tc 
write  (*,18)  fsam 
read  (*,*)  fsam 
write  (*,19)  nav 
read  (*,*)  nav 
nrep-  nav/100000 
ncount—  ifix  (1.0e7/fsam)  -  1 
call  adset  (ncount) 
write  (*,20) 
read  (*, '(a)')  filcal 
open  (3, file-f ileal, status-' old' ) 
read  (3,*)  ncals 
21    read  (3, ' (a)')  card 

if  (card.ne. '/*' )  goto  21 
avc  -  O.OdO 
do  25  k-1,32 
read  (3,*)  n,cal(k) 
25    avc  -  avc  +  cal(k) 
close  (3) 
avc  -  avc/32. 
do  27  k-1,32 

do  28  1-ncaI -ncals, ncal-1 
28    call(l,k)-  cal(k) /ncals 

do  27  1-0, ncal -ncals -1 
27    call(l,k)  -  O.OdO 
c 

c         start  big  loop  ; 
c 

nstept-  0 
nsp-  - 1 

write  (filnm(ilen+l: ilen+4) , ' (a) ' )  '.log' 
open  (4, file=filnm,status-'new' ) 
write  (4,123) 
do  300  11-1, ns 
nsd-  0 

peak  -  .false, 
peakh  -  . false . 
do  301  k-1,32 
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peakl(k)  -=  .false. 

301  peaklp(k)  -=  .true, 
do  302  k=l,31 
peak2(k)  -  .false. 

302  peak2p(k)  -  .true. 

write  (filnin(ilen+l:ilen+4)  ,23)  11+1000, dot 
write  (*,35)  filnm 
open  (3 .file-filnm, status- 'new' ) 
write  (3,*)  '#  of  averages  :  ' ,nav 
c        measure  f  and  q  in  transmission  : 

write  (1,*)  'output  01;A1B2A3' 

write  (1,7)  ampl(l:9) 

rewind  1 

ff  -  f 

df-  1.5*f/(parmt(2)*np) 

do  40  k-l,np/2 

if  (ll.eq.l)  then 
ff  -  ff  -  df 
else 

ff  =  ff  +  mult*df 
endif 
bb(k)  -  ff 

40  call  getp  (bb(k) , free ,aa(k) ) 

41  if  (ll.eq.l)  then 

ff  ■=  ff  -  df 
else 

ff  -  ff  +  mult*df 

endif 
bb(np/2+l)  -  ff 

call  getp  (bb(np/2+l),frec,aa(np/2+l)) 
if  (aa(np/2+l) .le.aa(np/2) .and.aa(np/2) . gt. 3 .e-7) 
1         goto  45 

do  42  k-l,np/2 
bb(k)-=bb(k+l) 

42  aa(k)-aa(k+l) 
goto  41 

45  do  48  k-np/2+2,np 

if  (ll.eq.l)  then 
ff  -  f f  -  df 
else 

f f  -  ff  +  mult*df 
endif 
bb(k)  -  ff 
48  call  getp  (bb(k) , free ,aa(k)) 

256  parmt  (1)  -  l.Od-30 

do  50  k-l,np 

parmt  (1)  -  dmaxl  (parmt(l) ,aa(k) ) 
50  if  (parmt(l)  .eq.  aa(k))  parmt(3)  -  bb(k) 

alamda  -  -1. 
55  chipr  -  chisq 

call  fit  (bb,aa,stdev,np, parmt, npar, chisq, alamda, funct) 

if  (abs(chisq-chipr)/ehisq  .gt.  0.0001)  goto  55 

dell  -  abs(fold  -  parmt (3)) 

if(ll.gt.l)  nstep-idnint  ( (del/dell)*nstep) 
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if  ((dell/del). It.. 001. and. 11. gt.l)  then 

write  (*,*)  'end  of  rod  travel  reached' 
goto  310 

endif 

nstep  -  min  (nstep,35) 

nstep  —  max  (nstep, 10) 

if  (11  .eq.  1)   fstart-  parmt(3) 
c         measure  f  and  q_hole  in  reflection  : 

write  (1,*)  'output  01;A2' 

write  (1,*)  'output  02;ap-20db' 

rewind  1 

fl  -  parmt(3)*(l.-20./parmt(2)) 

fr  -  parmt(3)*(l.+20./parmt(2)) 

call  getp  (fl, free, pi) 

call  getp  (fr,frec,pr) 

pin  -  0.5*  (pl+pr) 

call  getp  (parmt(3) , frec,pref ) 

refl  -=  1.  -  ((pin-pref)/pin) 

refl  -  10.*alogl0  (refl) 

write  (3,70)  parmt(3) ,parmt(2) 

write  (*,70)  parmt(3) ,parmt(2) 

write  (3,69)  refl 

write  (*,69)  refl 
c         check  cavity  temperature 

call  tem  (t) 

if  (t.gt.tc)  then 
write  (*,710) 
717  write  (*,720)  bell 

goto  717 
endif 

write  (3,74)  t 

write  (*,74)  t 
c       start  taking  data 

omega=  parmt(3) 

nlo-  0 

nsp-  nsp  +  1 

nrp-  mod(nsp,ncal) 
777  write  (*,84) 

nlo-  nlo  +  1 

if  (nlo  .gt.  nmax)  goto  778 

do  798  k-1,32 

798  spein(k)-0 

do  799  l-l,nrep 

call  spectrum  (speina, free ,parmt(3) ) 

do  799  k-1,32 

799  spein(k)-  spein(k)+speina(k) 
if  (nlo.eq.l)  then 

do  88  k-1,32 

88  skk(k)  -  dble(spein(k)) 
endif 

do  87  k-3,30 
87  spel(k)-  avc*(dble(spein(33-k))/cal(33-k)) 

do  89  k-3,29 

89  spe2(k)-  0.5*(  spel(k)+spel(k+l)  ) 
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aw  -  O.OdO 
do  827  k-3,30 
827  aw  -  aw  +  spel(k) 

aw  -  acvw/l  8  . 
do  90  k-3,30 

devl(k)   -   (spel(k) -aw)/(aw/dsqrt(dble(nav))) 
devs(nlo,k)    -  devl(k) 

if  (devl(k) .gt.5. )  peakh  -  .true, 
if  (devl(k) .ge.2. .and.peaklp(k))  then 
peak  -  . true . 
peakl(k)  -  .true. 

90  endif 
do  490  k-3,29 

dev2(k)-  sqrt(2.)*(spe2(k)-aw^)/(aw/dsqrt(dble(nav))) 
if  (dev2(k) .gt.5.)  peakh  -  .true. 
if  (dev2 (k). ge.2. .and. peak2p(k))  then 

peak  —  . true . 

peak2(k)  -  .true. 

490  endif 

write  (*,720)  bell 
write  (*,720)  bell 
call  plot  (devl ,omega-14. 8d-6 ,peakl ,peak) 
call  delay  (70) 
call  setmod  (3) 
if  (peak)  then 

write  (3,*)  'nlo=' ,nlo 

do  91  k-3,30 

if  (peakl(k))  then 

write  (*,720)  bell 

write  (*,92)  k,devl(k) ,omega-17 .d-6+(k-l)*l .093d-6 

write  (3,92)  k,devl(k) ,omega-17 .d-6+(k-l)*l .093d-6 

91  endif 

do  491  k-3,29 

if  (peak2(k))  then 

write  (*,720)  bell 

write  (*,492)  k,k+l,dev2(k) ,omega-16.5d-6+(k-l)*1.093d-6 

write  (3,492)  k,k+I ,dev2(k) ,omega-16 . 5d-6+(k-l)*l .093d-6 

491  endif 

peak-  . false . 
do  93  k-3,30 
peaklp(k)-  peakl(k) 

93  peakl(k)-  .false, 
do  94  k-3,29 
peak2p(k)-  peak2(k) 

94  peak2(k)-  .false. 
write  (*,97) 
goto  777 

endif 
778  if  (.not. peakh)  then 

do  888  k-1,32 
888  call(nrp,k)-  skk(k) 

else 

nsp-  nsp  -  1 

endif 
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avc"=  O.OdO 
do  892  k-1,32 
892  cal(k)-  O.OdO 

do  890  k-1,32 
do  891  l-0,ncal-l 
891  cal(k)-  cal(k)  +  call(l,k) 

890  avc=  avc  +  cal(k) 

avc-=  avc/32. 

write  (3, '(a)')  'calibration  file  :  '.filcal 
write  (3,*)  '#  of  spectra  :  ' ,nlo 
write  (3,103) 
do  895  k=3,30 
if  (nlo.eq.l)  then 

write  (3,104)  k,parmt(3) -17 .d-6+(k-l)*l . 093d-6 , 
1  devs(l,k) 

elseif  (nlo.eq.2)  then 

write  (3,105)  k,parmt(3) -17 .d-6+(k-l)*l .093d-6 , 
1  devs(l,k),devs(2,k) 

elseif  (nlo.eq.3)  then 

write  (3,106)  k,parnit(3) -17 .d-6+(k-l)*l . 093d-6 , 
1  devs(l,k) ,devs(2,k) ,devs(3,k) 

else 

write  (3,107)  k,parnit(3) -17  .d-6+(k-l)*l  .093d-6  , 
1  devs(l,k) ,devs(2,k) ,devs(3 ,k) ,devs(4,k) 

endif 
895  continue 

close  (3) 

if  (11  .eq.  ns)  fstop-  pannt(3) 
f  -  parnit(3)  -  mult*5*df 
fold  =  pannt(3) 
do  120  k=l,nstep 
if  (mult.eq.-l)  call  put  (1,1) 
if  (mult.eq.  1)  call  put  (0,1) 
120  call  delay  (1) 

nstept—  nstept  +  nstep 
write  (4,125)  ll,nlo,peakh 
300   continue 

310   write  (4,126)  f start, f stop, nstept 
rate—  l.e6*(fstop-fstart)/nstept 
write  (4,129)  rate 
close  (4) 

write  (filnin(ilen+l:  ilen+4)  ,  '  (a) ' )  '  .cal' 
open  (3  ,  file-f  ilnni,status-'new' ) 
write  (3,*)  ncal 
write  (3,103) 

write  (3,112)  (k,cal(k),  k-1,32) 
close  (3) 
stop 

7  format  ('output  02;', a) 

8  format  (f5.1) 

9  format  ('OEnter  LO  amplitude  in  transmission  (',f5.1,'  dB)  :  '$) 

10  format  ('OEnter  approx.  resonance  frequency  (  f  >  ftrue  )  :  '$) 

12  format  ('OEnter  approx.  value  of  Qloaded  :  '$) 

13  format  ('OEnter  #  of  stepper  motor  increments  :  '$) 
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69 

format 

70 

format 

72 

format 

74 

format 

80 

format 

82 

format 

84 

format 

92 

format 

492 


120 
130 
140 


format 


97 

format 

103 

format 

104 

format 

105 

format 

106 

format 

107 

format 

112 

format 

123 

format 

125 

format 

126 

format 

1 

129 

format 

710 

format 

720 

format 

end 

'$) 
'$) 


'$) 


'$) 
'$) 


dB') 

GHz',/,'  Q_loaded=  ',f7.0) 

K') 


f  5  .  2 , '  sigma' , 


12,'  with  ',f5.2,'  sigma', 


'OEnter  direction  of  rod  travel  (In  or  Out) 

'OEnter  filename  for  output  (no  extension) 

'OEnter  #  of  spectra  :  '$) 

'OEnter  max.  #  of  repeated  spectra  (',il,') 

'OEnter  critical  temperature  ( ' , f 3 . 1 , '  K)  : 

'OEnter  sampling  frequency  (',f7.0,'  Hz) 

'OEnter  #  of  averages  (',16,')  :  '$) 

'OEnter  filespec  for  filter  data  :  '$) 

i4,tl,al) 

'  spectrum  ' ,a) 

'  R_loss-  ' ,f6.2, 

'  F_res-  ' ,fll.9, 

'  Q_h/Q=.  ',f6.3) 

'  T_cavity=  ' , f 5 . 2 , 

flO.8) 

'output  02; ' ,a) 

/'  NOW  TAKING  DATA  '/) 

'  peak(l)  in  bin  #  ',12,'  with 

'   f-  ' , f 8 . 6 , '  GHz ' ) 
'  peak(2)  in  bin  #  ',12,'+' 
f-  ' , f 8 . 6 , '  GHz ' ) 
/'  DO  ANOTHER  SPECTRUM  '/) 
'/*') 

i2,lx,f9.7,lx,f6.2) 
i2,lx,f9.7,lx,f6.2,lx,f6.2) 
12  ,  Ix ,  f  9  .  7 ,  Ix ,  f  6  .  2  ,  Ix ,  f  6  .  2  ,  Ix ,  f  6  .  2 ) 
12 , Ix , f 9 . 7 , Ix , f 6 . 2 , Ix , f 6 . 2 , Ix , f 6 . 2 , ly  f 6  .  2 ) 
i2,3x,el3.6) 

'spectrum  averages  peakh') 
i6,5x,i2,5x,ll) 

'f(start)=  ',f8.6,'  GHz '/' f (stop)-   ' , f 8 . 6 , '  GHz'/ 
'nstept-  ' , 15) 

'tuning  rate=  ' , f 6 . 2 , '  kHz/step') 

///'  TEMPERATURE  TOO  HIGH  ,  MAGNET  QUENCH  POSSIBLE  !!!') 
'+',al) 


subroutine  getp  (f,frec,pow) 

real*8  f,pow 

character  frec*20,power*20 

write  (frec(3:12),120)  f 

write  (1,130)  frec(l:14) 

rewind  1 

call  delay  (3) 

write  (1,*)  'output  13;tr2' 

write  (1,*)  'enter  13' 

rewind  1 

read  (2)  power 

rewind  2 

read  (power(l : 11) , 140)  pow 
format  (flO.8) 
format  ('output  02;', a) 
format  (ell. 4) 
return 
end 
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subroutine  tem  (t) 

character*20  res 
write  (1,*)  'output  11;  tO?' 
write  (1,*)  'enter  11' 
rewind  1 
read  (2)  res 
rewind  2 

read  (res(l: 11) , 100)  r 
call  temp  (r,t) 
100   format  (ell. 3) 
return 
end 

subroutine  spectrum  (spein, free ,parm) 
c  program  to  start  DSP  program  and  upload  spectrum 
real*8  parm 
integer*4  spein(32) 
character  frec*20 
write  (frec(3:12) ,80)  parm+0.0107 
write  (1,82)  free 
write  (1,*)  'output  01;B1B3' 
write  (1,*)  'output  02;apl0.5db' 
rewind  1 
call  delay  (3) 
call  start 
call  tra  (64, spein) 
80      format  (flO.8) 
82      format  ('output  02;', a) 
return 
end 
subroutine  plot  (dev, fa,peakl ,peak) 
common  /windo/  ip(4) 
conraion  /cplot/  iw(4) 
common  /chsiz/  icwd,  ieht 
common  /tkl/  ixtkl(2) 

logical  xlog,ylog,peakl(32) ,peak,peak2(31) 
integer  xl , yb , xr , y t , cwd , cht , xtkl ( 2 ) 
real*8  dev(32),fa 

real*4  sca(2) ,scb(2) ,arr(28) ,bin(28) 
character  label*30 

equivalence  (iw(l) ,xl) , (iw(2) ,yb) , (iw(3) ,xr) , (iw(4) ,yt) 
data  sca(l)  /I./,  sca(2)  /I./,  scb(l)  /O./,  scb(2)  /O./ 
label(l:10)  -  'F  (kHz)  -  ' 
write  (label(ll:18),6)  fa 
6     format  (f8.6) 

label  (19:23)  -  '  GHz_' 
do  11  i  -  1,28 
arr(i)  -  sngl(dev(i+2)) 
11    bin(i)  -  (float(i)-l.)*1.093 
call  setmod  (18) 
iyaxis  -  yt  -  yb 
ixaxis  -  xr  -  xl 
xlo  -  0.0 
xhi  -  30. 
ylo  -  -10. 
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yhi  -  10. 

call  axis  (xl.yb, label , ixaxis , l,0,xlo,xhi , 30, 5) 

call  axis  (xl.yt,'   ', -ixaxis ,1,0, xlo.xhi , -30,5) 

call  axis  (xl ,yb, 'SIGMA_' , iyaxis ,0, l.ylo.yhi , -10, -2) 

call  axis  (xr,yb,'  ', -iyaxis,0,l,ylo, yhi, 10,2) 

xlog  -=  .false. 

ylog  -  . false . 

call  vtdata  (bin,xlog,xlo,xhi ,arr,ylog,ylo,yhi ,sca, scb,28 , 1, 12) 

if  (peak)  then 

call  cursor  (0,0) 

do  13  k-3,30 

13  if  (peakl(k))  write  (*,113)  k 

113  format  ('  PEAK  in  bin  #  ',12) 
do  14  k-3,29 

14  if  (peak2(k))  write  (*,114)  k,k+l 

114  format  ('  fat  PEAK  in  bin  #  ',12, '+',12) 
end  if 

call  cursor  (0,0) 

return 

end 

BLOCK  DATA  BLODAT 

C 

C  Block  data  subroutine  to  initialize  plotting  parameters  for 
C  program  PLOTRG. 

C--- - 

COMMON  /WINDO/  IP (4) 

COMMON  /CPLOT/  IW(4) 

COMMON  /CHSIZ/  CWD.CHT 

COMMON  /TKL/  XTKL(2) 

INTEGER  XTKL 

INTEGER  CWD,CHT 
C 

DATA  IP  70,0,639,479/ 

DATA  IW  780,40,620,400/ 

DATA  CWD  /8/,  CHT  /8/ 

DATA  XTKL  /10,5/ 

END 
************************************************************************ 


NAME  ADSET 
.MODEL  LARGE 
.CODE 

PUBLIC  ADSET 
ADSET  PROC 

push  bp  ;  establish  BP  as  frame  pointer 


program  to  set  sampling  rate  of  ADC. 
Callable  from  FORTRAN  {Call  ADSET(N) ) 
where  N  is  count  for  10  MHz  clock 
f_s  -  10  MHZ/  (N-1) 


mov  bp , sp 

push  es 

les  bx, [bp+6]  ;  get  count  from  caller 

mov  ax,es : [bx] 

mov  bx , ax 

mov  dx,0220h+0Bh  ;  define  Chimera  address 

mov  al,36h  ;  write  control  word  for  counter  0 

out  dx , al 

dec  dx 

dec  dx 


102 


pointing  to  0220h+8h 


write  count 


dec  dx 

mov  ax,bx 

out  dx.al 

mov  al , ah 

out  dx , al 

pop  es 

mov  sp.bp 

pop  bp 

ret  04h 
ADSET  ENDP 
END 
*********************************************************************** 


program  to  start  CHIMERA 
no  arguments 


enable  RI.O.RI.l  and  RESET 


.MODEL  LARGE 

.CODE 

PUBLIC  START 
START  PROC 

push  bp 

mov  bp , sp 

mov  dx,0220h 

inc  dx 

inc  dx 

mov  al , 1 

out  dx , al 

inc  dx 

mov  al, 10000000b 

out  dx , al 

dec  dx 

in  al ,  dx 

dec  dx 

in  al ,  dx 

inc  dx 

in  al ,  dx 

mov  al , 0 

out  dx.al 

in  al ,  dx 

inc  dx 

mov  al.lOOOOOOlb 

out  dx.al 

mov  sp , bp 

pop  bp 

ret 
START  ENDP 
END 
*********************************************************************** 


release  HOLD 


start  processor 


TRA 


NAME  TRA 
.MODEL  LARGE 
.CODE 

PUBLIC  TRA 
proc 

push  bp 

mov  bp , sp 

push  di 

push  si 

push  ds 


subroutine  to  transfer  data  from 
CHIMERA  to  host  .  To  be  called  by 
FORTRAN  program  {Call  TRA(N, ARRAY) ) 
where  N  is  #  of  words  and  ARRAY  is  of 
type  INTEGER 
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;  load  count  into  count  register  CX 

;load  destination  into  ES:DI 
; load  source  into  si 


; enable  transfer  from  Chimera 


check  HOLD  BIT 


enable  mem  space 


les  bx.dword  ptr[bp+10] 

mov  cx.es : [bx] 

les  di.dword  ptr[bp+6] 

mov  ax.OdOOOh 

mov  ds , ax 

mov  si,01800h*2 

mov  dx,0220h 

inc  dx 

inc  dx 
proc4:  in  al,dx 

mov  ah , 8 

and  al , ah 

cmp  al , 8 

je  proc4 
proc2 :  in  al , dx 
mov  ah , 1 

and  al , ah 

cmp  al , 1 

je  proc2 

mov  al, 11010000b 

out  dx , al 

push  ax 
procl:   in  al.dx 

and  al , 1 

jnz  procl 

pop  ax 

out  dx.al 

rep  movsw        ;  now  move  data 
mov  al, 11010001b 
out  dx,al 

pop  ds 

pop  si 

pop  di 

mov  sp.bp        ;  restore  BP 

pop  bp 

ret  08h 
TRA  ENDP 
end 

*  TMS320C25  assembly  program  for  32  point  realtime  spectrum  analysis. 

.sect  "init" 
RESET   B       INI      ;  branch  to  init  routine 
RIO     B       SERV 
*********************************************************** 

data  register 

Channel  0  control  register 

define  a  few  memory  Iocs 


data  page  0 

set  overflow  mode 

load  1  into  ACC 


into  host  memory 


SAMPLE  EQU 

0      ; 

CTRLO   EQU 

4      ; 

MEMO    EQU 

60h     ; 

MEMl    EQU 
*************v 

61h 
t-*********** 

.text 

INI   LDPK 

0      ; 

SOVM 

LACK 

1      ; 
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;   SET  IMR  (enable  INT  0) 

load  ARO  with  permanent  2 

load  ARl  with  address  of  averaged  sp 

load  AR2  with  address  of  first  buffer 

zero  accumulator 

point  to  AR2 

zero  1.  buffer 
point  to  ARl 


init.  aver,  spectrum 

;    start  A/D 

use  for  bigg  loop  control 
AR4  points  to  2.  buffer 
AR5  contains  count 
enable  interrupt 
discard  1.  sample 
AR4  points  to  2.  buffer 
AR5  contains  count 


use  for  big  loop  control 

page  6  and  1 .  buffer 

point  to  first  imag 

zero  routine 

do  32  point  power  spectrum 

AR2  points  to  new  spectrum 

average  result 

check  for  last  sample  ,  uses  mem  loc  0 

AR4  points  to  1.  buffer 

AR5  contains  count 

page  7  and  2 .  buffer 

point  to  first  imag 

zero 

do  spectrum 

AR2  points  to  new  spectrum 

do  average 

wait  for  last  sample  ,  uses  mem  loc  0 

2.  buffer 

count 

point  to  big  loop  counter 

branch 

point  to  bigg  loop  counter 

branch 


disable  interrupts 

transfer  aver,  spectrum  to  external  mem 


tell  host  I'm  finished 


************************************************************* 


SACL 

4 

LRLK 

0,2 

LRLK 

1 , 200h 

LRLK 

2 , 300h 

ZAC 

LARP 

2 

RPTK 

127 

SACL 

*+     ; 

LARP 

1       ; 

RPTK 

63 

SACL 

*+     ; 

SACL 

MEMO 

OUT 

^  ^  a.  a.  ^  a.  a,  ^  o.  o.  •4-4' •4"4' s 

MEMO , CTRLO 

-il' it- ■4' •^  •1"4' •4-4*  •!— 1- -4- 

LRLK 

j\«x  «3or«««3C« 

AR7,499 

LRLK 

AR4,380h 

LRLK 

AR5,1 

EINT 

CALL 

TEST 

LRLK 

AR4,380h 

LRLK 

AR5,64 

****** 

BIGGL  LRLK 

AR6,249 

BIGL  LDPK 

6 

LRLK 

AR2,301h 

CALL 

ZIMA 

CALL 

FFT64 

LRLK 

AR2,300h 

CALL 

AVER 

CALL 

TEST 

LRLK 

AR4 , 300h 

LRLK 

AR5,64 

LDPK 

7 

LRLK 

AR2,381h 

CALL 

ZIMA 

CALL 

FFT64 

LRLK 

2,380h 

CALL 

AVER 

CALL 

TEST 

LRLK 

AR4 . 380h 

LRLK 

AR5,64 

LARP 

6 

BANZ 

BIGL,*- 

LARP 

7 

BANZ 

BIGGL, *- 

***** 

DINT 

;  c 

CALL  OUTP 

;  t 

LDPK 

0 

LACK 

1 

SACL 

MEMO 

OUT 

MEMO , OCh 

IDLE  B 

IDLE 
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*******  interrupt  server  routine  **************************** 


SERV  DINT 

SST  MEMO 

LARP  4 

IN  *0+, SAMPLE, 5 

SBRK  1 

LDPK  0 

LST  MEMO 

EINT 

RET 
************************************************************** 


disable  interrupt 

save  status  reg 

point  to  buffer 

get  sample  and  store  in  buffer 

subtract  from  count  reg 

load  data  page  0 

restore  status  reg 

enable  interrupt 


TEST 


SAR 
LAC 
BGZ 
RET 


5,0 

0 

TEST 


store  (AR5)  in  0 

load  into  ACC 

wait  for  last  sample 


*************************************************************** 
ZIMA  LARP      2  ;  point  to  imags 

ZAC 

RPTK      63 

SACL      *0+  ;  zero  them 

RET 
**************************************************************^^^^^^ 

******   SPECTRUM  AVERAGE  ROUTINE  ************************** 


AVER 


AVLO 


l,0200h 
3,31 

;  point  to  av.  spectrum 
count 

1 
*+ 

*-,2 

point  to  average 
load  ACC 

*+,l 
*+ 

*+,0 
AVLO 

3 
*- 

1 

add  spectrum 
store 

(high  word  of  32  bit) 


LRLK 
LRLK 
LARP 
ZALH 
ADDS 
ADDS 
SACH 
SACL 
BANZ 
RET 
*****************************************************************^^ 

*********     OUTPUT  TO  EXTERNAL  MEM   **************************** 


OUTP 


LRLK  l,0200h 
LRLK  2,0201h 
LRLK  3,01800h 
LRLK  4,01801h 
LRLK  5,31 
1 


point  to  high  word  in  internal  mem 

low  word 

point  to  low  word  in  external  mem 

high  word 

count 


LARP 
OUS   LAC       *0+  ,0,4 

SACL      *0+,0,2 

LAC       *0+,0,3 

SACL      *0+,0,5 

BANZ      0US,*-,1 

RET 
*************** 

*****************************************************************^^ 
* 
*    Cooley-Tukey  Radix-4,  DIP  FFT  Program  -  64-point  straight-line. 


Three  radix-4  butterflies  -  implemented  with  macros. 
Complex  input  data  located  on  page  0  of  data  memory. 
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*  Results  are  left  in  data  RAM. 

*  Uses  13-bit  coefficients  from  MPYK  instructions. 

*  No  scaling  is  done  on  intermediate  values  in  the  program. 

*  Scaled  by  64 
* 

*  Data  memory  allocation. 
* 

XO  EQU  0          ;  X's  are  real  data  values 

YO  EQU  1           ;  Y's  are  imaginary  data  values 

XI  EQU  2 
Yl  EQU  3 
X2  EQU  4 
Y2  EQU  5 
X3  EQU  6 
Y3  EQU  7 
X4  EQU  8 
Y4  EQU  9 
X5  EQU  10 
Y5  EQU  11 
X6  EQU  12 
Y6  EQU  13 
X7  EQU  14 
Y7  EQU  15 
X8  EQU  16 
Y8  EQU  17 
X9  EQU  18 
Y9  EQU  19 
XIO  EQU  20 
YIO  EQU  21 

XII  "  EQU  22 
Yll  EQU  23 
X12  EQU  24 
Y12  EQU  25 
X13  EQU  26 
Y13  EQU  27 
X14  EQU  28 
Y14  EQU  29 
X15  EQU  30 
Y15  EQU  31 
X16  EQU  32 
Y16  EQU  33 
X17  EQU  34 
Y17  EQU  35 
X18  EQU  36 
Y18  EQU  37 
X19  EQU  38 
Y19  EQU  39 
X20  EQU  40 
Y20  EQU  41 
X21  EQU  42 
Y21  EQU  43 
X22  EQU  44 
Y22  EQU  45 
X23  EQU  46 
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Y23 

EQU  47 

X24 

EQU  48 

Y24 

EQU  49 

X25 

EQU  50 

Y25 

EQU  51 

X26 

EQU  52 

Y26 

EQU  53 

X27 

EQU  54 

Y27 

EQU  55 

X28 

EQU  56 

Y28 

EQU  57 

X29 

EQU  58 

Y29 

EQU  59 

X30 

EQU  60 

Y30 

EQU  61 

X31 

EQU  62 

Y31 

EQU  63 

X32 

EQU  64 

Y32 

EQU  65 

X33 

EQU  66 

Y33 

EQU  67 

X34 

EQU  68 

Y34 

EQU  69 

X35 

EQU  70 

Y35 

EQU  71 

X36 

EQU  72 

Y36 

EQU  73 

X37 

EQU  74 

Y37 

EQU  75 

X38 

EQU  76 

Y38 

EQU  77 

X39 

EQU  78 

Y39 

EQU  79 

X40 

EQU  80 

Y40 

EQU  81 

X41 

EQU  82 

Y41 

EQU  83 

X42 

EQU  84 

Y42 

EQU  85 

X43 

EQU  86 

Y43 

EQU  87 

X44 

EQU  88 

Y44 

EQU  89 

X45 

EQU  90 

Y45 

EQU  91 

X46 

EQU  92 

Y46 

EQU  93 

X47 

EQU  94 

Y47 

EQU  95 

X48 

EQU  96 

Y48 

EQU  97 

X49 

EQU  98 

Y49 

EQU  99 

X50 

EQU  IOC 
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Y50  EQU  101 

X51  EQU  102 

Y51  EQU  103 

X52  EQU  104 

Y52  EQU  105 

X53  EQU  106 

Y53  EQU  107 

X54  EQU  108 

Y54  EQU  109 

X55  EQU  110 

Y55  EQU  111 

X56  EQU  112 

Y56  EQU  113 

X57  EQU  114 

Y57  EQU  115 

X58  EQU  116 

Y58  EQU  117 

X59  EQU  118 

Y59  EQU  119 

X60  EQU  120 

Y60  EQU  121 

X61  EQU  122 

Y61  EQU  123 

X62  EQU  124 

Y62  EQU  125 

X63  EQU  126 

Y63  EQU  127 

Tl  EQU  7Eh        ;  Temporary  locations  addressed  ■=  126 

T2  EQU  7Fh        ;  by  auxiliary  registers.        -  127 

* 

*  Cosine  and  sine  values  (13-bit)  for  coefficients. 
* 

COO  EQU  4095 

COl  EQU  4076 

C02  EQU  4017 

C03  EQU  3919 

C04  EQU  3784 

COS  EQU  3612 

C06  EQU  3405 

C07  EQU  3166 

COS  EQU  2896 

C09  EQU  2598 

COlO  EQU  2275 

coil  EQU  1930 

C012  EQU  1567 

C013  EQU  1188 

C014  EQU  799 

C015  EQU  401 

COl 6  EQU  0 

C017  EQU  -401 

C018  EQU  -799 

C019  EQU  -1188 

C020  EQU  -1567 
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C021  EQU  -1930 

C022  EQU  -2275 

C023  EQU  -2598 

C024  EQU  -2896 

C025  EQU  -3166 

C026  EQU  -3405 

C027  EQU  -3612 

C028  EQU  -3784 

C029  EQU  -3919 

C030  EQU  -4017 

C031  EQU  -4076 

C032  EQU  -4095 

C033  EQU  -4076 

C034  EQU  -4017 

C035  EQU  -3919 

C036  EQU  -3784 

C037  EQU  -3612 

C038  EQU  -3405 

C039  EQU  -3166 

C040  EQU  -2896 

C041  EQU  -2598 

C042  EQU  -2275 

C043  EQU  -1930 

C044  EQU  -1567 

C045  EQU  -1188 

e046  EQU  -799 

C047  EQU  -401 

C048  EQU  0 

C049  EQU  401 

G050  EQU  799 

C051  EQU  1188 

C052  EQU  1567 

C053  EQU  1930 

C054  EQU  2275 

C055  EQU  2598 

C056  EQU  2896 

C057  EQU  3166 

C058  EQU  3405 

C059  EQU  3612 

C060  EQU  3784 

C061  EQU  3919 

C062  EQU  4017 

C063  EQU  4076 
* 

510  EQU  0 

511  EQU  401 

512  EQU  799 

513  EQU  1188 

514  EQU  1567 

515  EQU  1930 

516  EQU  2275 

517  EQU  2598 

518  EQU  2896 

519  EQU  3166 


no 


SIIO 

sill 

SI12 
SI13 
SI14 
SI15 
SI16 
SI17 
SI18 
SI19 
SI20 
SI21 
SI22 
SI23 
SI24 
SI25 
SI26 
SI27 
SI28 
SI29 
SI30 
SI31 
SI32 
SI33 
SI34 
SI35 
SI36 
SI37 
SI38 
SI39 
SI40 
SI41 
SI42 
SI43 
SI44 
SI45 
SI46 
SI47 
SI48 
SI49 
SI50 
SI51 
SI52 
SI53 
SI54 
SI55 
SI56 
S157 
SI58 
SI59 
SI60 
SI61 
SI62 
SI63 


EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU  - 

EQU  - 

EQU  - 

EQU  - 

EQU  - 

EQU  - 

EQU  - 

EQU  - 

EQU  - 

EQU  - 

EQU  - 

EQU  - 

EQU 


3405 

3612 

3784 

3919 

4017 

4076 

4095 

4076 

4017 

3919 

3784 

3612 

3405 

3166 

2896 

2598 

2275 

1930 

1567 

1188 

799 

401 

0 

-401 

-799 

-1188 

-1567 

-1930 

-2275 

-2598 

-2896 

-3166 

-3405 

-3612 

-3784 

-3919 

-4017 

-4076 

-4095 

-4076 

-4017 

-3919 

-3784 

-3612 

-3405 

-3166 

-2896 

-2598 

-2275 

-1930 

1567 

1188 

799 

401 


Ill 


*********************************************************************** 

* 

*  MACROS  TO  PRODUCE  STRAIGHT-LINE  64-POINT  COMPLEX  FFT. 
* 

*********************************************************************** 
* 

*  ZERO  for  case  theta  -  0. 
* 

*  X's  and  Y's  are  input  and  output  locations  for  butterfly. 


XII  -  XII  +  XI3 
XI3  =  XII  -  XI3 


ZERO     $MACR 

D  XI, X 

LAC 

:XI1 

ADD 

:XI3 

SACL 

:XI1 

SUB 

:XI3 

:,1 

SACL 

:XI3 

LAC 

:XI: 

SUB 

:XI2: 

SACL  ■ 

fr.O.ARl 

ADD 

:XI2: ,1 

ADD 

:XI1: 

SACL 

:XI: 

SUB 

XII  :,1 

SACL 
* 

XI2: 

LAC 

YIl 

ADD 

YI3 

SACL 

YIl 

SUB 

YI3 

,1 

SACL 

YI3 

LAC 

YI: 

SUB 

YI2: 

SACL  ' 

t,0,AR2 

ADD 

YI2:,1 

ADD 

YIl: 

SACL 

YI: 

SUB 

YI1:,1 

SACL 

YI2: 

LAC  -i 

Ir 

ADD  : 

YI3 

SACL 

XII 

SUB 

YI3 

,1 

SACL  ■* 

t.O.ARl 

LAC  ^ 

^0,AR2 

ADD 

XI 3 

SACL 

YI3 

SUB 

XI3. 

,1 

SACL  : 

YIl: 

LAC  ■* 

r 

SACL  : 

XI3: 

Tl 

-  XI  - 

XI2 

Rl 

(ACC) 

=  XI  + 

XI 2 

XI 

=  Rl  + 

XII 

XI2 

-  Rl 

-  XII 

YIl  -  YIl  +  YI3 
YI3  -  YIl  -  YI3 

T2  -  YI  -  YI2 

Rl  (ACC)  -  YI  +  YI2 

YI  -  Rl  +  YIl 

YI2  -  Rl  -  YIl 

XII  -  Tl  +  YI3 
Tl  -  Tl  -  YI3 

YI3  -  T2  +  XI3 
YIl  -  T2  -  XI3 
XI3  -  Tl 


LAC 


:XI: ,15 
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SACH 

XI: 

LAC 

XII 

,15 

SACH 

XII 

LAC 

XI2 

,15 

SACH 

XI 2 

LAC 

XI 3 

,15 

SACH 

XI 3 

LAC 

YI: 

15 

SACH 

YI: 

LAC 

YIl 

,15 

SACH 

YIl 

LAC 

YI2 

,15 

SACH 

YI2 

LAC 

YI3 

,15 

SACH 

YI3 

$ENDM 
************************************************************* 


* 

* 

* 

* 

* 

NORM 

* 


NORM  -  standard  radix-4  butterfly. 

X's  and  Y's  are  locations  of  data  input  and  output. 
lA's  specify  twiddle  factor  locations. 


$MACRO  XI ,XI1 ,XI2 ,XI3 , YI , YIl , YI2 , YI3 , lAl , IA2 , IA3 


LAC 
ADD 
SACL 
■SUB 
SACL 
LAC 
SUB 


XII 
XI 3 
XII 
XI 3 
XI 3 
XI: 
XI 2 


,1 


SACL  *,0,AR1 


ADD 

ADD 

SACL 

SUB 

SACL 


:XI2; 

:XI1: 

:XI: 

:XI1; 

:XI2: 


XII  -  XII  +  XI3 


XI 3 

-  XII 

-  X 

Tl 

-  XI  - 

XI 2 

Rl 

(ACC)  - 

=  XI 

+  XI2 


XI 


Rl  +  XII 


XI2  -  Rl  -  XII 


LAC 

ADD 

SACL 

SUB 

SACL 

LAC 

SUB 


YIl 
YI3 
YIl 
YI3 
YI3 
YI: 
YI2 


SACL  *,0,AR2 


ADD 

ADD 

SACL 

SUB 

SACL 


:YI2: 

:YI1: 

:YI: 

:YI1: 

:YI2: 


YIl  -  YIl  +  YI3 


YI3  -  YIl  -  YI3 


T2  -  YI  -  YI2 

Rl  (ACC)  =  YI  +  YI2 


YI  -  Rl  +  YIl 


YI2  -  Rl  -  YIl 
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LAC 

* 

ADD 

:YI3: 

SACL 

:XI1: 

SUB 

:YI3:,1 

SACL 

*,0,AR1 

LAC 

* 

ADD 

:XI3: 

SACL 

* 

SUB 

:XI3: ,1 

SACL 

:YI1: 

LT 

:YI1: 

MPYK 

C0:IA1: 

PAC 

LT 

:XI1: 

MPYK 

SI:IA1: 

SPAC 

LT 

:YI1: 

SACH 

:YI1: ,4 

MPYK 

SI:IA1: 

PAC 

LT 

:XI1: 

MPYK 

COrlAl: 

LTA 

:YI2: 

SACH 

:XI1: ,4 

MPYK 

C0:IA2: 

PAC 

LT 

:XI2: 

MPYK 

SI:IA2: 

SPAC 

LT 

:Y12: 

SACH 

:YI2: ,4 

MPYK 

SI:IA2: 

PAC 

LT 

:XI2: 

MPYK 

CO : IA2 : 

LTA 

*,AR2 

SACH 

:XI2:,4 

MPYK 

CO : IA3 : 

PAC 

LT 

*,AR1 

MPYK 

SI:IA3: 

SPAC 

SACH 

:YI3:,4 

MPYK 

CO : IA3 : 

PAC 

LT 

*,AR2 

MPYK 

A  PAP 

SI:IA3: 

ArAU 

SACH 

:XI3:,4 

XII  -  Tl  +  YI3 
Tl  -  Tl  -  YI3 

T2  =  T2  +  XI3 
YIl  -  T2  -  XI3 


YIl  =  YI1*C01  -  XI1*SI1 


XII  -  YI1*SI1  +  XI1*C01 


YI2  =  YI2*C02  -  XI2*SI2 


XI2  -  YI2*SI2  +  XI2*C02 


;  YI3  -  T2*C03  -  T1*SI3 


;  XI3  -  T1*C03  +  T2*SI3 


LAC 


:XI: ,15 
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SACH 

XI: 

LAC 

XII 

,15 

SACH 

XII 

LAC 

XI2 

,15 

SACH 

XI 2 

LAC 

XI 3 

,15 

SACH 

XI 3 

LAC 

YI: 

15 

SACH 

YI: 

LAC 

YIl 

,15 

SACH 

YIl 

LAC 

YI2 

,15 

SACH 

YI2 

LAC 

YI3 

,15 

SACH 

YI3 

$ENDM 
*********************************************************************** 
* 

*  SPEC  for  case  theta  -  pi/4 
* 

*  X's  and  Y's  are  data  input  and  output  locations. 
* 

SPEC 
* 


$MACRO  XI ,XI1 ,XI2 ,XI3 , YI , YIl , YI2 , YI3 


LAC 

ADD 

SACL 

SUB 

SACL 

LAC 

ADD 

SACL 

SUB 

SACL 

LAC 

ADD 

SACL 

SUB 

SACL 

LAC 

ADD 

SACL 

SUB 

SACL 

LAC 

SUB 

SACL  *,0,AR1 

ADD   :XI:,1 

SACL  :XI: 


XI: 

XI2: 

XI: 

XI2: 

XI2: 

YI: 

YI2: 

YI: 

YI2: 

YI2: 

XII 
XI 3 
XII 
XI 3 
XI 3 

YIl 
YI3 
YIl 
YI3 
YI3 

:XI1; 
:XI: 


;  XI  -  XI  +  XI2 
;  XI2  -  XI  -  XI2 

;  YI  =  YI  +  YI2 
;  YI2  -  YI  -  YI2 

;  XII  -  XII  +  XI3 
;  XI3  -  XII  -  XI3 

;  YIl  -  YIl  +  YI3 
;  YI3  -  YIl  -  YI3 

;  Tl  =  XII  -  XI 
;  XI  -  XII  +  XI 
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LAC 

YI: 

ADD 

YIl 

SACL 

YI: 

SUB 

YIl 

,1 

SACL  * 

LAC 

XI 2 

ADD 

YI3 

SACL 

XII 

SUB 

YI3 

,1 

SACL 

YI3 

LAC  *,0,AR2 

SACL 

XI2 

LAC 

YI2 

SUB 

XI 3 

SACL 

YIl 

ADD 

XI3 

,1 

SACL 

XI 3 

LAC  i 

t 

SACL 

YI2 

LT 

YIl 

MPYK  ( 

:o8 

PAC 

LT 
MPYK  C 

XII 

:o8 

SPAC 

SACH 

YIl 

.4 

APAC 

LTA 

YI3 

SACH 

XII 

,4 

ZAC 

MPYK  C 

;o8 

SPAC 

LT    : 

XI3: 

MPYK  C 

)08 

SPAC 

SACH  : 

YI3: 

,4 

APAC 

ArAO 

SACH  : 

XI3: 

,4 

LA 
SA 

C 
lCH 

XI: 

XI: 

15 

LA 

SA 

c 

CH 

XII 
XII 

,15 

LA 
SA 

C 
CH 

XI 2 
XI 2 

,15 

LA 
SA 

C 
CH 

XI 3 
XI 3 

,15 

LA 

SA 

C 
CH 

YI: 
YI: 

15 

LA 
SA 

C 
CH 

YIl 
YIl 

,15 

YI  -  YI  +  YIl 
T2  -  YI  -  YIl 

XII  -  XI2  +  YI3 
YI3  -  XI2  -  YI3 
XI2  =  T2 

YIl  -  YI2  -  XI3 
XI3  =  YI2  +  XI3 
YI2  =  Tl 


YIl  -  (YI1-XI1)*C08 


XII  -  (YI1+XI1)*C08 


;  YI3 


(YI3+XI3)*C08 


;  XI3  -  <XI3  -  YI3)*C08 
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LAC 

YI2 

,15 

SACH 

YI2 

LAC 

YI3 

,15 

* 

* 

SACH 

YI3 

- 

$ENDM 

BITRVI 
* 

$MACRO     PR , QR 

* 

ZALH 

PR 

store  Re(l)  in  ACCH 

ADDS 

QR 

store  Re(2)  in  ACCL 

SACL 

PR 

store  ACCL 

SACH 

QR 

store  ACCH 

ZALH 

PR 

+1 

ADDS 

QR 

+1 

SACL 

PR 

+1 

SACH 

QR 

+1 

$ENDM 

********************************************************************** 
* 

* 

MAIN  routii 

le  to  call 

above  macros  with  appropriate  parameters. 

************************************************************************ 


FFT64 


ZAC 

LRLK  1,1 
LARP  1 
RPTK  63 
SACL  *0+ 
LARP  2 
LRLK  AR2,T1 
LRLK  AR1,T2 


zero  ACC 
point  to  imags 

store  0  in  imags 

AR2  points  to  Temp 
ARl  points  to  Temp 


*  PASS 
* 


ZERO  X0,X16,X32,X48,Y0,Y16,Y32,Y48 

NORM  XI ,X17 ,X33 ,X49 , Yl , Y17 , Y33 , Y49 ,1,2,3 

NORM  X2,X18,X34,X50,Y2,Y18,Y34,Y50,2,4,6 

NORM  X3,X19,X35,X51,Y3,Y19,Y35,Y51,3,6,9 

NORM  X4,X20,X36,X52,Y4,Y20,Y36,Y52,4,8,12 

NORM  X5 ,X21 ,X37 ,X53 , Y5 , Y21 , Y37 , Y53 , 5 , 10, 15 

NORM  X6,X22,X38,X54,Y6,Y22,Y38,Y54,6,12,18 

NORM  X7,X23,X39,X55,Y7,Y23,Y39,Y55,7,14,21 

SPEC  X8,X24,X40,X56,Y8,Y24,Y40,Y56 

NORM  X9 ,X25 ,X41 ,X57 , Y9 , Y25 , Y41 , Y57 , 9 , 18 , 27 

NORM  X10,X26,X42,X58,Y10,Y26,Y42,Y58,10,20,30 

NORM  Xll ,X27 ,X43 ,X59 , Yll , Y27 , Y43 , Y59 , 11 , 22 , 33 

NORM  X12,X28,X44,X60,Y12,Y28,Y44,Y60,12,24,36 

NORM  X13,X29,X45,X61,Y13,Y29,Y45,Y61,13,26,39 

NORM  X14,X30,X46,X62,Y14,Y30,Y46,Y62,14,28,42 
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NORM   X15,X31,X47,X63,Y15,Y31,Y47,Y63,15,30,45 


*  PASS  2 


ZERO  X0,X4,X8,X12,Y0,Y4,Y8,Y12 

ZERO  X16,X20,X24,X28,Y16,Y20,Y24,Y28 

ZERO  X32 .X36 ,X40,X44 , Y32 , Y36 , Y40 ,Y44 

ZERO  X48,X52,X56,X60,Y48.Y52,Y56,Y60 

* 

NORM  XI ,X5 ,X9 ,X13 , Yl , Y5 , Y9 , Y13 ,4 , 8 , 12 

NORM  X17 ,X21 ,X25 ,X29 , Y17 . Y21 , Y25 , Y29 ,4, 8 , 12 

NORM  X33,X37.X41,X45,Y33,Y37,Y41,Y45,4,8,12 

NORM  X49 ,K53 ,X57 ,X61 , Y49 , Y53 , Y57 , Y61 ,4, 8 , 12 

* 

SPEC  X2,X6,X10,X14,Y2,Y6,Y10,Y14 

SPEC  X18 ,X22 ,X26 ,X30 , Y18 , Y22 , Y26 , Y30 

SPEC  X34,X38,X42,X46,Y34,Y38,Y42,Y46 

SPEC  X50,X54,X58,X62,Y50,Y54,Y58,Y62 

* 

NORM  X3 ,X7 ,X11 ,X15 , Y3 ,Y7 . Yll , Y15 , 12 , 24 , 36 

NORM  X19 ,X23 ,X27 ,X31 , Y19 , Y23 ,Y27 , Y31 , 12 , 24, 36 

NORM  X35 ,X39 ,X43 ,X47 , Y35 , Y39 , Y43 , Y47 , 12 , 24 , 36 

NORM  X51,X55,X59,X63,Y51,Y55,Y59,Y63,12,24,36 

* 

*  PASS  3 
* 

ZERO  X0,X1,X2,X3,Y0,Y1,Y2,Y3 

ZERO  X4,X5,X6,X7,Y4,Y5,Y6,Y7 

ZERO  X8,X9,X10,X11,Y8,Y9,Y10,Y11 

ZERO  X12,X13,X14,X15,Y12,Y13,Y14,Y15 

ZERO  X16 ,X17 ,X18 ,X19 , Y16 , Y17 ,Y18 , Y19 

ZERO  X20,X21,X22,X23,Y20,Y21,Y22,Y23 

ZERO  X24,X25,X26,X27,Y24,Y25,Y26,Y27 

ZERO  X28 ,X29 ,X30,X31 , Y28 ,Y29 , Y30, Y31 

ZERO  X32 ,X33 ,X34,X35 , Y32 ,Y33 , Y34, Y35 

ZERO  X36 ,X37 ,X38 ,X39 , Y36 , Y37 ,Y38 , Y39 

ZERO  X40,X41,X42,X43,Y40,Y41,Y42,Y43 

ZERO  X44  ,X45  ,X46  ,X47  ,  Y44 ,  Y45  ,  Y46  ,  Y47 

ZERO  X48 ,X49 ,X50,X51 , Y48 , Y49 , Y50 , Y51 

ZERO  X52,X53,X54,X55,Y52,Y53,Y54,Y55 

ZERO  X56 ,X57 ,X58 ,X59 , Y56 ,Y57 ,Y58,Y59 

ZERO  X60,X61,X62,X63,Y60,Y61,Y62,Y63 
* 

*  UNSCRAMBLE  OUTPUT 
* 

BITRVI  2,32 
BITRVI  4,64 
BITRVI  6,96 
BITRVI  10,40 
BITRVI  12,72 
BITRVI  14,104 
BITRVI  18,48 
BITRVI  20,80 
BITRVI  22,112 
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BITRVI 

26,56 

BITRVI 

28,88 

BITRVI 

30,120 

BITRVI 

36,66 

BITRVI 

38,98 

BITRVI 

44,74 

BITRVI 

46,106 

BITRVI 

52,82 

BITRVI 

54,114 

BITRVI 

60,90 

BITRVI 

62,122 

BITRVI 

70,100 

BITRVI 

78,108 

BITRVI 

86,116 

BITRVI 

94,124 

*** 

*  COMPUTE  POWER  SPECTRUM 

*** 

SQU     $MACRO  ADR.ADI.ADS 

SQRA     :ADR 

SQRA     :ADI 

SACH     :ADS 

ZAC 

* 

$ENDM 

******************* 

ZAC 

LT 

0 

* 

MPY 

0 

SQU 

1,2,0 

SQU 

3,4,1 

SQU 

5,6,2 

SQU 

7,8,3 

SQU 

9,10,4 

SQU 

11,12,5 

SQU 

13,14,6 

SQU 

15,16,7 

SQU 

17,18,8 

SQU 

19,20,9 

SQU 

21,22,10 

SQU 

23,24,11 

SQU 

25,26,12 

SQU 

27,28,13 

SQU 

29,30,14 

SQU 

31,32,15 

SQU 

33,34,16 

SQU 

35,36,17 

SQU 

37,38,18 

SQU 

39,40,19 

SQU 

41,42,20 

SQU 

43,44,21 

SQU 

45,46,22 

ACC  ->  ACC  +  P,  [ADR]* [ADR]  ->  P 


load  real  0  into  T 
multiply  and  store  in  P 
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SQU  47,48,23 
SQU  49,50,24 
SQU  51,52,25 
SQU  53,54,26 
SQU  55,56,27 
SQU  57,58,28 
SQU  59,60,29 
SQU  61,62,30 
SQRA  63 
APAC 
SACH      31, 

****************** 
RET 

********** 
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